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A new Monte-Carlo based uncertainty analysis is introduced to quantita- 
tively determine the predictive ability of QCD sum rules. A comprehensive 
analysis of ground state p-meson and nucleon spectral properties is performed. 
^ I Many of the findings contradict the conventional wisdom of both practitioners 

Ch ' and skeptics alike. Associations between the phenomenological fit parame- 

ters are particularly interesting as they reveal how the sum rules resolve the 
^ . spectral properties. The use of derivative sum rules for the determination 

^ \ of /9-meson spectral properties is shown to be a very unfavorable approach. 

Most prior nucleon sum rule analyses are based on a sum rule which is found 
to be invalid; the results are suspect, and should be reevaluated. The "loffe 
formula", argued by many to qualitatively encapsulate a description of the 
nucleon mass in terms of the chiral symmetry breaking order parameter i^q) 
is misleading at best. QCD Sum Rules are found to be self-consistent without 
contributions from direct instantons. This implies that instanton effects are 
adequately accounted for in the nonperturbative vacuum condensates. This 
in-depth examination of QCD sum rule self consistency paints a favorable 
picture for further quantitative refinements of the QCD sum rule approach. 
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I. INTRODUCTION 



A. Prologue 

The QCD Sum Rule (QCD-SR) approach to QCD continues to be a highly active field. 
The infiuence of the field is refiected in over 1500 references to the seminal paper of Shifman, 
Vainshtein and Zakharov (SVZ) [|[], with over 380 of these references following the beginning 
of 1992 While the approach has been applied to a variety of hadronic observables in 
both vacuum and finite density nuclear matter, a comprehensive review of the systematic 
errors associated with the approach has not been considered. 

One of the key assumptions of the approach is the use of the so called "continuum 
model" to remove excited state contaminations from the hadron correlator under investi- 
gation. Recently, the continuum model for the nucleon was tested by using correlation 
functions calculated via lattice regularized QCD The results suggest that the contin- 

uum model adequately removes the excited state contaminations and allows the isolation of 
the ground state. The success of the continuum model warrants a more careful investigation 
of other aspects of the QCD-SR approach so that uncertainties in the predictions may be 
reliably and quantitatively determined. 

The focus of this paper is to establish a rigorous procedure for extracting quantities of 
phenomenological interest from QCD Sum Rules. The main conclusion of this investigation 
is that QCD Sum Rules work and they are predictive when the analysis is done rigorously. 
The ground state masses of the p-meson and the nucleon are used to introduce the concepts. 
However, the procedures introduced here are more generally applicable to quantities of 
current experimental interest. The goal is to shed light on the stability of QCD Sum Rule 
analyses and thus the accuracy to which observables may be reliably determined. 

Implementation of the procedure introduced here does not lead to slight adjustments 
in the extracted parameters, but rather leads to qualitatively different conclusions. It will 
become apparent that most existing QCD-SR investigations of nucleon properties are in 
fact unreliable. Previous calculations are based on sum rules in which neither the opera- 
tor product expansion (OPE) nor the phenomenological description are under acceptable 
control. 

The QCD-SR method is presently the best fundamentally based approach for investi- 
gating the properties of hadrons in nuclear matter, as the lattice approach is challenged by 
a number of formidable obstacles 0. Issues surrounding the finite width of the rho- meson 
make reliable lattice calculations of p — a; mixing very difficult. These are two of the many 
topics of current interest to nuclear physicists where the best fundamentally based approach 
is that of QCD-SRs. Moreover, the approach has minimal model dependence. Hence it is 
imperative to probe the predictive ability of the approach and the ideas presented here are 
currently being applied to these topics mentioned above 

B. Systematic Uncertainties 

The most significant source of uncertainty in the QCD Sum Rule approach is an imprecise 
knowledge of the vacuum condensates appearing in the operator product expansion (OPE). 
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This is particularly problematic for the higher dimension operators, where one usually in- 
vokes the vacuum saturation hypothesis (factorization) to replace higher-dimension-operator 
vacuum-expectation values by products of lower- dimension operators. A list of key sources 
of uncertainty in the QCD-SR approach should include: 

• the unfactorized condensate values, 

• factorization of higher dimensional operators, 

• small but neglected corrections, 

• the truncation of the operator product expansion (OPE), 

• the selection of the regime for matching the QCD and phenomenological sides of the 
sum rules, 

• uncertainties associated with the summation of perturbation theory in large orders, 
and 

• the possibility of significant direct instanton contributions. 

With the exception of direct instanton contributions, the effects of these uncertainties will 
be estimated via a Monte-Carlo error analysis. 

Gaussian distributions for the condensate values are generated via Monte Carlo. These 
distributions are selected to reflect the spread of values assumed in previously published 
QCD-SR analyses, and the uncertainties such as operator factorization listed above. These 
distributions provide a distribution for the OPE and thus uncertainty estimates for the OPE 
which will be used in a fit. In turn, the OPE distribution provides distributions for the 
phenomenological fit parameters, thus establishing the predictive ability of the QCD-SR 
approach. 

C. Conventional Analyses 

In the field of QCD inspired models, the focus is often on whether or not the model is able 
to encompass the known experimental data, as opposed to assessing the predictive potential 
of the method. A review of the QCD-SR literature leaves no doubt as to whether or not 
the QCD-SR method is able to encompass the established data. The poorly determined 
condensate values and arbitrary aspect of the matching Borel regime provides tremendous 
freedom in devising ways to achieve the desired result. In fact, it is difficult to find QCD- 
SR-based predictions of hadronic observables which fail to agree with experiment. This 
is somewhat surprising, since the physics responsible for some observables may be poorly 
represented in a truncated OPE commonly limited to distances of less than 0.3 fm. 

Most QCD-SR calculations of nucleon properties in the literature fall short in the analysis 
stage where the phenomenological parameters are determined by matching the QCD and 
phenomenological sides of the sum rules. Typical analysis short-comings include: 

• selecting a single value for the Borel parameter which gives "nice" results. A Borel 
regime should be selected in order to evaluate the stability and reliability of the results. 
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• a selection of the Borel regime without careful regard to OPE convergence^] nor main- 
taining ground state dominance of the phenomenological side of the sum rules. 

• the fixing of search parameters (such as the continuum threshold) to preferred val- 
ues. It will become apparent that this introduces a strong bias to the remaining fit 
parameters which may not reflect the properties of QCD. 

• claiming an accuracy for QCD-SR predictions without supporting calculations. Occa- 
sionally a "stability analysis" 0JT^ is considered, in which fit parameters are monitored 
as a single condensate value is varied. However, such analyses explore a relatively small 
corner of the condensate parameter space. 

In general, former procedures adopted for matching the two sides of the sum rules lack the 
level of rigor which will be presented here. 

This general lack of rigor can be traced back to the original paper of SVZ |^. In their 
conclusions, they comment "we prefer not to deepen into the debris of computer calcula- 
tions, sophisticated fit programs, arranging error bars here and there. . ." While this may be 
appropriate for a seminal paper on the QCD-SR approach, continuation of this philosophy 
has lead to strong negative criticism of the QCD-SR field. It is not uncommon to hear 
remarks such as "QCD Sum Rules are useful only if you know the answer." or "You can get 
anything you want from QCD Sum Rules." and "10% here, 10% there. . . Pretty soon you're 
talking about real numbers!" These remarks are unfortunate, as a rigorous and proficuous 
analysis is possible. 



D. Outline 

In this paper, a method for the quantitative determination of phenomenological quanti- 
ties with uncertainties will be presented. Since this paper is directed to both QCD-Sum-Rule 
practitioners and skeptics alike, we begin by briefly outlining the QCD-SR approach in Sec- 
tion A more pedagogical review of the approach may be found in the recent review 
article of Ref. |rT[. This section also introduces the "loffe formula" IT^, argued by many to 



qualitatively encapsulate a description of the nucleon mass in terms of the chiral symmetry 
breaking order parameter (gg). 

Uncertainties associated with the evaluation of the OPE will be addressed via a new 
Monte-Carlo based procedure. Estimates for these uncertainties are presented in Section 
[m| . Section |^ details the Monte-Carlo uncertainty analysis and discusses the optimization 
algorithm and cautious convergence criteria used in the following. 

Section presents an analysis of p-meson sum rules. Here the criteria for the selection of 
the regime for matching the QCD and phenomenological sides of the sum rules is reviewed. 
Reasonable alternatives are considered via Monte Carlo. 



^Here and in the following, "convergence" of the OPE simply means that the highest dimen- 
sion terms considered in the OPE, with their Wilson coefficients calculated to leading order in 
perturbation theory, are small relative to the leading terms of the OPE. 
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These considerations lead to the selection of the optimal interpolating field for nucleon 
sum rules. New insights into unconventional nucleon interpolators obtained from the lattice 
QCD investigation of Ref. |T^ play a significant role here. Section |VJ reviews the selection 
of the optimal nucleon interpolating field for correlators obtained from spin-1/2 interpolat- 
ing fields. The effects of choosing non-optimal interpolators is demonstrated. Section |V11| 
explores additional sum rules obtained from the overlap of the generalized spin-1/2 and a 
spin-3/2 interpolator. 

In Section [VII1| we explore the correlations between condensate values, fit parameters 
and sum rule consistency. These correlations serve to identify the roles of the leading terms 
of the OPE in hadronic physics. The validity of the "loffe formula" is evaluated here, and 
these results may be of interest to those modeling the QCD vacuum. The contingency table 
analysis reveals that the intimate relationship between the quark condensate and the nucleon 
mass suggested in the "loffe formula" is invalid. 

The necessity of direct instanton contributions to the QCD Sum Rules is determined 
in Section We demonstrate that there is no evidence indicating that direct instanton 
contributions are required to maintain sum rule consistency. 

The implications of the results of this analysis on coordinate space correlators are dis- 
cussed in Section ^ Finally, Section ^ summarizes the findings of this in-depth investiga- 
tion. 



II. QCD SUM RULE FORMALISM 

The underlying principle of field theoretic approaches to hadron phenomenology is the 
Ansatz of duality. That is, it is possible to simultaneously describe a hadron as quarks 
propagating in the QCD vacuum, and as a phenomenological field with the appropriate 
quantum numbers. 

Hadron masses are extracted from an analysis of the two-point function 

U{p) =z Jd'x e'P-' {Q I T{xix) x(0)} | Q) . (2.1) 

The interpolating field x is typically constructed from quark field operators combined to 
give the quantum numbers of the hadron under investigation. The p-meson interpolator 
related to the decay constant is 

= ^ (^7m u - dj^ d) . (2.2) 

To maintain maximal overlap with the ground state relative to excited states, only interpo- 
lators without derivatives are considered. For the nucleon we begin by considering the most 
general spin-1/2 nucleon interpolator 

Xo = Xi + f^X2 , (2.3) 

where 

X^{x) = e"^" (u^''ix)C-f5d\x)) u'ix) , (2.4a) 
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and 



(2.4b) 



Xi{x) is the interpolator typically used in lattice QCD analyses. X2{.x) vanishes in the 
nonrelativistic limit. However in a theory with light relativistic current quarks, there is no 
reason to exclude such an interpolating field a priori |]14|. With the use of the Fierz relations, 
the combination of the above two interpolating fields for (3 = —1 may be written 



= 2 (X2 - Xi) , 



(2.5a) 
(2.5b) 



giving the proton interpolating field advocated by loffe and often found in QCD-SR 
calculations. 



A. Phenomenology 

At the phenomenological level one proceeds by inserting a complete set of eigenstates 
with the quantum numbers of the interpolator. The ability of the interpolator to annihilate 
a positive parity baryon to the QCD vacuum is described by the parameter A in 

(f2 I xo(0) I s) = AoiM(p, s) , (2.6) 

for state i of momentum p and spin s. u{p, s) is a Dirac spinor. Negative parity states 
require a 75 on the right-hand side of ( |2.6| ). For the nucleon, the two point function has the 
form 

n(p) = EAa. ^,_^._^^. , (2.7) 

where +/— corresponds to positive /negative parity baryons. While there is some suppression 
of excited state contributions to the two-point function, there is little hope of isolating ground 
state properties from such a function. Hence one uses the celebrated Borel Transform 

B[f{p')]= lim ^{-P'r^'(-f^] fip'), (2.8) 

n— »oo ^ ' 

-p'^/n=M'^ 

at each Dirac-7 structure to obtain exponential suppression of excited states. The parameter 
M is commonly referred to as the Borel mass. The Borel transform is easily applied to a 
dispersion relation for H(p) at each Dirac-7 structure 

^/ 9x 1 /""^ , IniH(s) 

H(p^) = - / ds V + subtractions , (2.9) 

vr JO s — p^ 

and the polynomial subtraction terms are eliminated by the Borel transform. At the struc- 
ture 7 ■ p one has 
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n^.p(M) = J p^.p{s) e-'/^^' ds , (2.10a) 

where we have introduced the spectral density p{s) = Imn(s)/7r, 

p^.p{s) = \l 6{s - Ml) + i,.,{s) . (2.10b) 

At the structure 1 

ni(M) = j pi(s) e""/^^' ds , (2.11a) 

where 

pi{s) = (±M^) Xl 6{s - Ml) + Us) ■ (2.11b) 



Here, ^7.p(s) and ii{s) account for excited state contributions, and +/— in (|2.11b|) corre- 
sponds to positive/negative parity states. 
For the vector meson one has 

nv(M) = / pv{s) e-'l^^ ds , (2.12a) 



Pvis) = f^Sis-M'^)+^yis), (2.12b) 
at the structure g^i, — V^iVvIv^ ■ 

B. QCD 

At the quark level, one exploits the operator product expansion (OPE) to describe the 
short distance behavior of the two-point function 

T{x(x), x(0)} =Y.Cn{x,p) a(0,/i), (2.13a) 

n 

= Co{x) I 

+ Ci{x) TUq 

+ Cs{x) qq 

+ C4.i(x) G^.G'^^^ + C,,2ix) m^qq 

+ C5.i(x) ga^.y gG"^'^ + C,,2{x) rnqGl^G^^' 

+ G6.i(a;) qVqqVq^ C^.^ix) m^qa qG'^^'' ^ C^,^{x) f a^cGl^Gl^Gl^ 
+ G,,^{x) qq Gl^G^^^ + G7.2(x) mqqVq qVq + G,,^{x) f a,cGl,G\^G\^ 
+ G8.i(x) qVq qTa^^^qG'^^" + Gs.2{x) m^qq G^G'^''" 
+ GsAx) Gl^G^^'G^yG^'P^ , (2.13b) 
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where /i is the normahzation point at which the coefficient functions and the operators are 
defined. Here we have exphcitly included all operators up to dimension eight, to leading 
order in the quark mass rrig, having the quantum numbers of the vacuum. The first digit of 
the subscript of the Wilson coefficients C{x) gives the energy dimension of the operator. F 
may take any of the 16 independent Dirac-7 matrices. 

Upon Fourier transforming to momentum space, ( |2.13|) becomes an expansion in l/Q"^, 
valid for large momentum transfers. Thus one may use perturbation theory to calculate the 
Wilson coefficients. 

The contribution of the gluon operator term 4.1 in the expansion for the proton is 
suppressed due a factor of (27r)^ which appears in the denominator of the Wilson coeffi- 
cient arising from integration over the quark-loop momentum which introduces the factor 
(i^g/(27r)^. Terms involving more than two gluon field strength tensors and products of 
nig and C^^^G""^^ are estimated to be small and are typically neglected. 

The standard treatment of the OPE proceeds via inserting the explicit forms of the 
interpolating fields into the two-point function of ( p.l|) and contracting out pairs of time- 
ordered quark-field operators which are the fully interacting quark propagators of QCD. 
Wick's theorem for the time-ordered product of two fermion fields provides 

(fi|T{g»,g''(0)} 

= (0 I T[q\x),q\Q)} \ 0) + |: g'^(x),g^(0) :| fi) . (2.14) 

Taylor expansion of the normal ordered piece leads to 

:g^(x),g^(0): = : gS(0), g^^(O) : + : 9,g^(0),g^^(0) : 

+ y^x'':d,d,ql{Q),%{Q): +■■■. (2.15) 

The well known trick is to select the coordinate gauge A^ = Q such that 

M^) = ~Ig,M x""-^ (9aG^.(0)) x^x-" + ■■■ , (2.16) 



and the partial derivatives of ( |2.15|) may be replaced by covariant derivatives 



d^^V^ = d^ + igA^. (2.17) 



Generally, one works with ( |2.16| ) to leading order. The following relations allow an easy 



determination of the covariant derivatives 

t'^V^C = (2.18a) 

[V^VI =z(7eG^% (2.18b) 

V'g^ = -^^7ca-Gg». (2.18c) 
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In summary, the quark correlators employed in calculating the OPE under the usual as- 
sumption of vacuum saturation of the intermediate states of composite operators include 

{n\T{q^{x),q\0)}\Q) = 

7 ■ X S""^ - 



m„ 



227r2a;2 22 3 

. 2 



mg{qq) 7 ■ x 5"'' + 7^{q9cCT ■ Gq)& 



ab 



2^3 



26 3 



2732 
i 



mg{qgc(T ■ Gq) 7 ■ x5 



ab 



X 



21033 



+ 



2^7r2x2 
1 



A 



{9cGl,){^-xa-^ + a-^^-xr-f 



In 



'-x^A^' 



+ 27 



EM 



I „ r^n \ nl3 ^ab 



(2.19) 



{n\T{q\x)g,Gl^i'{Q)]\n) 

1 , ^ . X 
"2 



263 

i 



{qgc<y ■ Gq)(y^p-'''^ 



28 ^^q{(19c(T ■ Gq) (7 ■xaap + Oo^p 1 ' 



X 



21032 



{glG'){qq)a^,^. 



(2.20) 



Finally the momentum space correlator is Borel transformed. Here one encounters the 
factorial suppression of higher-dimension operators as indicated in 



B 



1 



1 



(A; - 1)! (M2) 



2\fc-l ■ 



(2.21) 



It is worth noting that the factorial suppression does not really set in until one reaches 
the term For the p meson this term corresponds to an operator of dimension eight, 

and for the nucleon the dimensions are twelve and thirteen for the two nucleon sum rules 
obtained from spin-1/2 interpolating fields. Since the OPE for the p-meson is traditionally 
truncated at dimension eight f^] , the correlator should be more reliable than for the nucleon 
sum rules which are truncated at dimension nine. 



C. Power Corrections from the Large-Order Behavior of Perturbation Theory 



It is well known that the standard treatment of the OPE encounters difficulties at large 
orders in the perturbative expansion of the Wilson coefficients. Perturbative expansions are 
divergent at large orders and are asymptotic at best [TB|. In the standard treatment, these 
divergences limit the accuracy to which the perturbative coefficients may be determined. 

Renormalons [|1^], a particular set of perturbative graphs, provide a simple illustration 



of the factorial growth in the contributions of perturbative graphs at large orders |[T8|-p2 



The factorial growth in large-order contributions associated with soft infrared (IR) virtual 
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momenta is problematic, as the series is not Borel summable. As such, the Wilson coeffi- 
cients, by themselves, are ill defined in the standard treatment. The uncertainty resulting 
from the restriction to finite orders of the asymptotic expansion is commonly referred to 
as the IR renormalon ambiguity. Of course the use of perturbative propagators to describe 
soft virtual momenta is incorrect. Instead, one must utilize fully dressed nonperturbative 
propagators. Hence, it is not surprising one encounters difficulties in the IR regime, as one 
is simply encountering the Landau pole of the perturbative coupling constant. 

A solution to the IR renormalon ambiguity problem was provided by the ITEP group 
through the introduction of the normalization/separation scale fi of (|2.13a|). The OPE 



becomes a separation of scales, with virtual momenta lying above /i represented in the 
Wilson coefficients, and momenta below /i in the vacuum expectation values (VEVs) of 
higher-dimension operators beginning at dimension four. The introduction of the scale n 
eliminates the IR renormalon ambiguity of the perturbative expansion. However, in princi- 
ple, both the Wilson coefficients and the VEVs of operators contain perturbative as well as 
nonperturbative physics. 

In principle, the OPE defined in ( p.l3a| ) is free of IR renormalon ambiguities. However, 
in practice it is difficult to implement the formal procedure of truncating virtual momenta 
below fi when calculating higher order corrections to the Wilson coefficients. While these 
contributions are small for leading order corrections, their inclusion is a source of error which 
can be accommodated in the Monte-Carlo uncertainty analysis. 

In the ultraviolet (UV) virtual momenta regime, the renormalon series is sign alternating 
and Borel summable. Summation of the UV series can give rise to a power correction 
proportional to Aqqj-j/Q^ which closely resembles the squared quark mass contribution to 
the OPE 1^ . For the vector correlator, the UV renormalon gives rise to an effective mass of 
roughly (2/7r)^/^ Aqcd/3, which is much larger than the usual current quark mass the order 
of 5 MeV [^. Hence there is some concern that important power corrections associated 
with UV renormalons has been overlooked in the standard QCD-SR treatment. Moreover, it 
is the power corrections that provide the link to nonperturbative phenomena. To the extent 
that the power correction is scheme dependent and not reliably known, we will refer to these 
quark-mass like power corrections as UV renormalon uncertainties. 

Additional evidence for 1/g^-type power corrections independent of renormalons has 
been obtained through a consideration of constraints imposed by asymptotic freedom and 



analyticity on the large-order behavior of perturbation theory ||2^ , |25[| . In this case the 
perturbative series is not Borel summable. Given the absence of a local gauge-invariant field- 
operator product of dimension two, some might argue that the coefficient of the singularity 
in the Borel plane governing the large-order behavior of perturbative series must have a 
zero. However, perturbative estimates of the coefficient do not provide evidence of such a 
p5| and the issue is unsettled. 



zero 



Hence the issue of possible dimension- two power corrections is of current interest. The 
techniques presented in this manuscript might be used to provide some phenomenological 



^For a particularly lucid discussion of why the standard treatment of the OPE is phenomenolog- 
ically successful, see Ref. [20|. 
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insight into the debate. However, such considerations would take us too far afield and 
will be deferred to a subsequent analysis [26|. On the other hand, it is interesting to look 



for discrepancies which might be rooted in such power corrections, and this is done in the 
following analysis. 



D. loffe Formula 



To summarize the points of this section, we present the traditional QCD Sum Rules for 
the nucleon [|l5l obtained from the consideration of the interpolating field Xsr/2 of (|2.5| ) to 



operator dimension 8. 



' ^ 8 (27r)' 32 {2nf n ' 6 

+ ^{q9(^-Gq) {qq)^ 



= A^e-^^-/^^' + --- . (2.22a) 

At 1 : --^-^{qq)M^ + —{qq){^G^) 
4 (27r)^ 288^^^^ ^TT ' 

= A^M;ve-^^/^'' + -- - . (2.22b) 

The dots allow for excited state contributions.^] 

The loffe Formula for the nucleon mass is obtained from (|2.22|) by keeping the leading 



terms of the OPEs and assuming dominance of the ground state on the right-hand side of 
the equations, 

l^M«^A^e-*^^/^% (2.23a) 

-\ ^(gg) ^ A^M^ve-*^^/^^^ . (2.23b) 

The loffe formula follows from the ratio of these two equations, 

M^^-lM^, ,2.24a, 
= 0.90 GeV for M = 1 GeV. 

If one prefers to eliminate the Borel mass by selecting M = Mm, 



^The Wilson coefficient of the dimension 7 operator in ( 2.22b ) has had an elusive history. Early 
calculations |1^,27,2|] appear to have calculated a subset of the diagrams contributing to the 
dimension 7 operator However a more persistent error lies in a factor of 3 correction in the 
term [x^/(2^° 3^)] {qq) {glC^) Oapy^-^jl of ( |2.20D from that used in previous analyses. 
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Mn-{-2 {271 f (qq))'^' , (2.24b) 
= 0.97 GeV. 

Here the quark condensate is estimated from the partially conserved axial current (PCAC) 
relation and pion decay to give 

(qq) = ^ , (2.25) 

= (0.225 GeV)^ , 

which corresponds to = .140 GeV, J't^ = 0.093 GeV and the average light quark mass is 
7 MeV. The qualitative features of chiral symmetry breaking coupled with the quantitative 
success of these relations has lead to enormous faith being placed into the validity of the 
loffe formula. Most practitioners expect this relation to be valid, at least in a qualitative 
manner of proportionality. 



E. Excited State Contributions 

The operator product expansion provides knowledge of the two-point function ( |2.1j ) in 
the near perturbative regime. As such the correlator represents contributions from both 
the ground state hadron and excited states. To maintain some predictive ability in the 
QCD-SR approach, a "continuum model" for excited state contributions to the correlator is 
introduced. 

Guided by the principle of duality, the leading terms of the OPE surviving in the limit 
M — > CX3 are used to model the form of the spectral density accounting for excited state 
contributions. Consider, as an example, the identity operator for the nucleon interpolator 
XsR of ( p^) . Exphcitly including the excited state term of ( p.lO|) in equation (|2.22a| ) leads 
to 

= / ^^.p{s) e""/*^ ds , (2.26) 



2 (2 7r)4 

for large Borel mass. The Laplace transform readily provides 

The continuum model contribution is obtained by introducing a threshold, Sq, at which 
strength in the excited states becomes significant. Such a model is motivated well by the 
experimental cross section for e"*" e~ — > hadrons. In addition, this approach accommodates 
the finite widths and the multi-particle continuum of QCD. Hence the excited state contri- 
butions associated with the identity operator are 

e-/- 4 ^ M«e--'M' . ^ . l) . (2.28) 
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Continuum model contributions associated with and terms are calculated in a similar 
manner. Usually the continuum contributions are placed on the OPE side of the sum rules. 
Hence, terms surviving in the limit M — oo have the following factors associated with 
them: 




(2.29a) 

(2.29b) 
(2.29c) 

The contributions of this model relative to the ground state, whose properties one is really 
trying to determine, are not small. They are typically 10 to 50% 0. The validity of this 
somewhat crude model is relied upon to cleanly remove the excited state contaminations. 

In the lattice QCD investigation of the continuum model formulated in Euclidean space 
0, the short time regime of point-to-point lattice correlation functions is described well by 
the QCD-Sum-Rule-inspired continuum model. There, the Laplace transform of the spectral 
density appears to be sufficient to render any structure in the spectral density insignificant 
in the short Euclidean time regime of point-to-point correlators. Similar conclusions are 
expected to hold for the Borel transformed sum rules under investigation here. 

It should be noted that the spectral density of the physical world is different than that 
encountered on the lattice. For example, in the nucleon channel one expects strength in the 

m 



whereas in the lattice 



correlator above the ground state to start at the vrp threshold 
calculations such contributions are altered by the heavier quark mass and the quenched 
approximation. While the lattice results are encouraging, one may find that in some cases 
the continuum model is not sufficiently detailed when the continuum model contributions 
are large. An example of this will be presented in the discussion surrounding the p-meson 
sum rules, where discrepancies may be due to an over simplified continuum model. 



III. UNCERTAINTY ESTIMATES 

The origin of the uncertainty estimate for the OPE lies in the uncertainties assigned 
to the QCD parameters appearing in the truncated OPE. In this section we assess the 
uncertainties on these QCD parameters. 

In assigning these uncertainties one would like to select values conservatively enough such 
that the QCD-SR approach can test the validity of our present understanding of QCD. QCD- 
SR predictions that fail to agree with experiment at the 2a level should indicate a possibly 
interesting discrepancy. At the same time one would also like to select uncertainties in accord 
with the generally accepted values in the literature. In addition, the error assessments need 
to refiect uncertainties associated with the factorization of higher dimensional operators and 
small but never-the-less neglected corrections. These considerations have been balanced 
in the following. While some may argue that some values are known better than that 
indicated here, others will no doubt argue that the errors are underestimated. 

In any event, one will learn for the first time how the uncertainties in the QCD parameters 
are mapped into uncertainties in the phenomenological fit parameters. In the lattice QCD 
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investigation of Ref. it was noted that interplay between the pole position and continuum 
threshold lead to rather large uncertainties in the fit parameters. It will be interesting to 
see how reasonable estimates of the QCD uncertainties are revealed in the fit parameters. 

For the quark condensate, (gg) we note that there are two values commonly found in 
QCD-SR analyses. PCAC and pion decay considerations coupled with quark masses of 
rriu + rrid = 11 MeV lead to a value of —(0.250)'^ GeV^. QCD-SR considerations of qa 
and octet baryon magnetic moments ppj-p^ prefer a smaller magnitude at —(0.225)^ GeV^. 
Expecting the value to lie some where in this regime, we take the average of these values with 
an uncertainty of half the difference at an 80% confidence level. Introducing the standard 
QCD-SR notation we define 



i27i f{qq) = 0.52 ± .05 GeV^ 



(3.1) 



This value and uncertainty is also in accord with a 
recent quark mass estimates |53 . 



0.47 ± 0.05 GeV^ obtained from more 



Early estimates of the gluon condensate (a^/TrG^) from charmonium sum rules H] place 
the value at 0.012 ± 0.005 GeV^, which is generally referred to as the "standard value" 



b = {27ry{— G^) = 0.47 ± 0.20 GeV^ . 

TX 



(3.2) 



However, a number of more recent investigations place the gluon condensate at much larger 



values p4|-[36[]. For example, a finite-energy sum rule (FESR) analysis of the p-meson channel 
places the gluon condensate at 2 to 5 times the standard value |3^. In this analysis the 
spectral density is taken from a fit of the isospin-one channel of e"*" e~ scattering data. First 



results from the ALEPH and CLEO II experiments l3Jj33] also suggest a larger value at 
{as/nG'^) = 0.025 ± 0.010 GeV^. A recent Laplace based sum rule analyses suggests a 
value of {os/tiG'^) ~ 0.015 ±0.005 GeV^, more in accord with the standard value. However, 
we are somewhat skeptical of the refinements of this result as the analysis is based on an 
approach which is criticized in the following.^ The most reliable sum rules of Sec. VII also 
prefer a larger gluon condensate in order to maintain consistency between the sum rules. 
There it is found that the nucleon mass may be reproduced with a value of {ag/n G"^) ^ 0.03 
GeV^ which agrees with experimental estimate. Since the gluon condensate plays a critical 
role in the nucleon sum rules as part of a factorized dimension seven operator we adopt an 
uncertainty of 50% and take 



b = {2n f{— G^) = 1.2 ± 0.6 GeV^ 



7r 



(3.3) 



The mixed condensate is parameterized as {q gc<^ ■ G q) = —Tn^{qq). Early analyses 
of baryons as well as more recent analyses of heavy-light quark systems place 



For example, while the language used to describe the sum rules in Ref. [38| is somewhat different 
from that used here, the nature of the sum rules are the same. The concerns surrounding OPE 
convergence, continuum models and the use of derivative sum rules discussed in this paper still 
apply. 
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mQ = 0.80 GeV^, while an analysis of N3/2 QCD-SR sum rules demand tuq to be smaller 
0] at ttiq = 0.65 GeV^. Hence we take the average and half the difference as the uncertainty. 

ml = 0.72 ± 0.08 GeV^ . (3.4) 

Since the parameter ttiq is selected via Monte Carlo, the value of the mixed condensate is 
obtained by multiplying ml by the central value of the quark condensate, as opposed to a 
randomly selected value of the quark condensate. This is in contrast to factorized operators 
where both condensates appearing in the operator are selected via Monte Carlo. 

Relatively little is known about the magnitude of the dimension-six, four-quark operators. 
Early arguments placed the values of these condensates within 10% of the factorized values 
^ . However other analyses claimed significant violation of factorization [Il0| , p4| - |36|j38|j41|j42 



for these operators in both nucleon and vector meson sum rules. Parameterizing the four 
quark-operators as k {qq)'^, estimates of factorization violation place k = 2 or more. For the 
nucleon we consider 

K7v = 2±l, and 1.0<K7v<4.0. (3.5) 

If one adopts the standard value for the gluon condensate, then Kp ~ is sufficient to 
reproduce the p-meson mass. However, selection of the more recent and reliable estimates 
of 6 ~ 1.2 GeV^ demands ~ 6 to reproduce the p-meson mass. This value is in accord 



with the observations of Refs. p^j35[] , however the uncertainty on this parameter is large. 
Hence we adoptj^ 

Kp = 6±2, and 1.0 < Kp < 10.0 . (3.6) 

Figure |l| displays histograms drawn from 1000 QCD parameter sets for four of the QCD 
parameters discussed to this point. 

Variation of the QCD scale parameter Aq^j-, has little effect on the results. However, we 
restrict the values to the conventional range adopted in QCD-SRs 

Ag)^j3 = 0.15 ± 0.04 GeV, and 0.10 GeV < Ag^^ < 0-20 GeV . (3.7) 

(3) 

Variation of Aq^j-, above 200 MeV is only sensible if the separation scale p of the OPE, taken 
to be 500 MeV, is increased to maintain plausible convergence of the power corrections of 
the OPE. A change in p redefines the condensate values. A systematic determination of the 
new condensate values is beyond the scope of this investigation and we choose to simply 
restrict the value as above. 

(3) 

It is perhaps interesting to note that if Aq^j^ is very large, say near the upper limit of 
the current world average |^^, the OPE separation scale would necessarily become large. 



Consequently, the Borel mass would be restricted to larger values, and access to ground state 



^With Kp ~ 6 as opposed to 2, even a small suppression of the chirally even four-quark condensate 
at finite density can give rise to substantial suppression of the p-meson mass and width. Further 
investigation of these issues is warranted. 
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FIG. 1. Histograms for QCD parameter distributions drawn from a sample of 1000 QCD 
parameter sets. Included are (a) the quark condensate parameter of ( |3.1| ), (b) the gluon condensate 
parameter of ( |3.3D , (c) the mixed condensate parameter mg and (d) the vacuum saturation violation 
parameter Kp. 



properties may not be possible. The correlator may be dominated by strength from excited 
states, similar to the present situation for the pseudoscalar pion correlator. Fortunately, 
the lower bound on the current world average for Aq^j-, bodes well for present QCD sum 

rule analyses. However, if the current world average for Aqqj-, persists, a comprehensive 
reformulation of the approach should be pursued. 

The p-meson sum rule requires knowledge of Os/vr at the scale of 1 GeV^. Since this 
scale is rather low for two-loop perturbation theory, we adopt a 10% uncertainty 



aJl GeV^l 



0.117±0.012. 



TT 



(3.8) 



Finally, the p-meson sum rule requires knowledge of the average light quark mass rrtq. 
Since it appears multiplied by the quark condensate, we will use the relation of ( p.25|) and 
equate 
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mq (qq) = - -ml fl . (3.9) 

Monte-Carlo techniques might also be used to probe OPE truncation errors by randomly 
selecting the size of the next term to follow the truncated series. The problem lies in 
determining a sensible and realistic estimate without actually calculating the term. Instead, 
we will emphasize the importance of OPE convergence and select the Borel analysis regime 
such that truncated terms are very unlikely to be important. 



IV. QCD-SR ANALYSIS 

A. Monte-Carlo Uncertainty Analysis 

Gaussian distributions are generated using established algorithms for the conversion of 
uniform distributions on the interval (0,1) |Q. A particularly fast algorithm used in this 
investigation is as follows. If Ui and U2 are uniform on (0,1), then construct fi = 2mi — 1 
and f2 = 2^2 — 1. Calculate = f ^ + v^- If > 1 start over. Otherwise, calculate 

/-2 1nr2y/^ ^ /-2 1nr2y/^ 
^i = vi\ — I ) and Z2 = V2\ — I . (4.1) 

zi and Z2 are independent and normally distributed with mean and variance 1. 

A set of Gaussianly-distributed randomly-selected condensate values is generated, from 
which an OPE in Borel space, 11'-'^^ (M), is constructed. By repeating this procedure, an 
uncertainty for the OPE may be determined. Selecting ub evenly distributed points in the 
Borel parameter space {Mj,j = 1, 2, . . . ns), the standard deviation in the OPE at the j'th 
Borel mass is given by 



'^OPE 



1 2 

(M,) = Y: n«P^(M,) - noPE(M,) , (4.2) 

nn — i - I L J 



nc 

Tie 

where i denotes the i'th set of nc QCD parameter sets and 



1 nc 

nopE(M,) = — Enr^(M,). (4.3) 

i=i 

In practice, nc = 500 was used in estimating the OPE uncertainty. However, nc = 100 was 
sufficient to get a stable uncertainty estimate. 

The uncertainties in the OPE are not uniform throughout the Borel regime. They are 
larger at the lower nonperturbative end of the Borel region where uncertainties in the higher- 
dimensional vacuum condensates dominate. Hence, it is crucial that the appropriate weight 
is used in the calculation of x^- Having determined the OPE uncertainty, a measure is 
easily constructed. For the OPE obtained from the fc'th set of QCD parameters, the per 
degree of freedom is 

Xl 1 [HO™(M,)-HSi^(M,;A,,m,,^,)]' 

o (4.4) 



Ndf riB-Upf^^ o-QPEiMj] 
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where Up is the number of phenomenological search parameters, and n^^(Mj; A^, m^, w^) 
denotes the phenomenological Spectral Representation of the QCD-SR. In practice, ub = 51 
points were used along the Borel axis, is minimized by adjusting the phenomenological 
parameters including the pole residue Aq, the hadron mass m, and the continuum threshold 
w. 

Distributions for the phenomenological parameters are obtained by minimizing for 
many QCD parameter sets. Provided the resulting distributions are Gaussian, an estimate 
of the uncertainty in the phenomenological results such as the hadron mass is obtained from 

I N N 

= ^ ~ ^^"^ ' J^H^k- (4.5) 

k=l k=l 

In the event that the resultant distribution is not Gaussian, we will report the median 
and asymmetric standard deviations from the median. We note that ([4. 5]) is the standard 



deviation of the distribution and not the standard error obtained by further dividing by 
A^, the number of QCD parameter sets in the sample. The standard deviation is roughly 
independent of the number of QCD parameter sets used in the analysis, and directly reflects 
the input parameter uncertainties. In the following, we generally select N = 1000. While 
100 QCD parameter configurations are sufficient for obtaining reliable uncertainty estimates, 
we also wish to explore correlations among the QCD parameters and the phenomenological 
fit parameters. The extra sets aid in resolving more subtle correlations. 

Some care must be taken in the interpretation of the x^/^df- While we are free to 
choose any number of points ub along the Borel axis for fitting the two sides of the sum 
rule, it is important to recognize that the actual number of degrees of freedom is determined 
by the number of terms in the OPE. For the x^/^df to have true statistical meaning, a 
correlated x'^ calculation is required. 

XI = T.T. [^T'^m - nSR(M,; A., m., wu) x 

[nOP^(M,0 - nSi^(M,v; Afc, mfc, wu)] , (4.6) 
where Cjj, is the inverse covariance matrix. The covariance matrix may be estimated by 

1 ■^c 

c,,, = ^^^--^ 5: [nr^(M,) - nopE(M,)] [n?p^(M,0 - nopE(M,,)] (4.7) 

Unfortunately, for most of the fits considered here the covariance matrix is ill-conditioned. 
Inverting the covariance matrix via singular value decomposition leads to well known patho- 



logical problems |^ such as best fit lines which do not pass through the data points etc. 
While one could experiment with different approaches to stepping around this problem, we 
prefer to use the robust method of determining best fit parameters by minimizing the uncor- 
related x^ of (Eil)- Even if some correlation among the data is suspected, it is still acceptable 
to use the uncorrelated x^ fit provided the errors on the parameters are estimated using a 
Monte-Carlo based procedure like that discussed here rather than from the dependence of 
X^ on the parameters pl,E5[. 
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When simultaneously fitting multiple sum rules which may have some fit parameters in 
common, ( [1.4|) is simply modified to 



DF nsR Ub - Up 



i=l j=l 



(4.8) 



where i runs over the considered sum rules. (In this example, m is a common parameter 
to all nsR sum rules.) The subscript i in Mij indicates that each sum rule has a unique set 
of hb Borel masses Mj spanning the region of interest in Borel space. The subscript k in 
( |4.4|) , denoting a particular set of QCD parameters, has been suppressed. In the following 
discussion it will become apparent that it is advantageous to weight the contributions in the 
sum over the various sum rules by the size of their valid Borel regimes. If is the size of 
the valid Borel regime for the i'th sum rule, then (|4.8|) is modified to 



X' 



1 



Ndf nsR Ub - n„ J2i '^'^ 



"SR. "S 

EE 

B i=l j=l 



V, 



B 



npPE(M,,)-nP(M,,;A„m,m 



t^OPEi(Mi) 



(4.9) 



B. Search Algorithm 

The optimization method utilized in this investigation is an updated version of a 



direction-set routine by Powell originally published in Ref. |^^. This algorithm finds the 
minimum value of a function by iterative variation of the function parameters which need 
not be independent. Modifications of the original routine improve the selection of conju- 
gate directions in the search and avoid the possible generation of linearly dependent search 



directions [47 



Powell's method is among the best methods for finding minima in a multidimensional 



parameter space when derivatives are not readily available P8| , ^ . Part of the reason for the 
robust nature of Powell's routine is the cautious criterion for ultimate convergence. Once a 
minimum is found, the optimal vector, a, in the parameter space is displaced by ten times 
the requested accuracy. The search proceeds from this new point until a minimum is found 
at vector h. The minimum on the line joining a and h is found at c. Convergence is reached 
if the components of the vectors (a — c) and {h — c) are all less than 10% of the required 
accuracy. Otherwise the direction (a — c) is utilized in a new search. 

The displacement of the first minimum by ten times the requested accuracy plays a 
crucial role in removing the sensitivity of the ultimate minimum from the initial parameter 
estimates. While there is sufficient flexibility in the algorithm parameters to search in the 
locality of the initial parameter estimates, the opposite criteria has been utilized in this 
investigation. One of the main goals of this investigation is to eliminate the possibility of 
fine tuning the initial parameter estimates which otherwise can bias the results. Instead, we 
are interested in the predictions of the QCD sum rules. 
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V. p-MESON SUM RULES 



The fundamental QCD-SR for the p meson follows from the consideration of ( ^.11 ), ( ^.21 ), 
( p.l2| ) and ( p.29|) . We consider the correlator, lisiM"^), where the continuum model is 



subtracted from the phenomenological side and placed on the OPE side. 



C2_ C3 



(m - l)r(M2 



The Wilson coefficients are given by |jT| 



1 ( <^s\ 

Co = ^ 1 + — , ci = 0, 
C2 = mg (gg) + (— G^^ , 

Z4 TT 

112 



C3 = -—T^a,Kp{qq) . (5.2) 

The phenomenology of the QCD-SR is described by the vector-meson pole of interest 
plus the continuum model accounting for the contributions of all excited states. By working 
in a region where the pole dominates the phenomenology, one can minimize sensitivity to the 
model and have assurance that it is the spectral parameters of the ground state of interest 
that are being determined by matching the sum rules. In practice, these considerations 
effectively set an upper limit in the Borel parameter space, beyond which the model for 
excited states dominates the phenomenological side. 

At the same time, the truncated OPE must be sufficiently convergent as to accurately 
describe the true OPE. Since the OPE is an expansion in the inverse squared Borel mass, 
this consideration sets a lower limit in Borel parameter space, beyond which higher order 
terms not present in the truncated OPE are significant and important. Monitoring OPE 
convergence is absolutely crucial to recovering nonperturbative phenomena in the sum-rule 
approach, as it is the lower end of the Borel region where the nonperturbative information 
of the OPE is most significant. As we shall see, this information must also be accurate. 

In short, one should not expect to extract information on the ground state spectral 
properties unless the ground state dominates the contributions on the phenomenological 
side and the OPE is sufficiently convergent. In this investigation, we will analyze each 
individual sum rule with regard to the above criteria. A sum rule with an upper limit in 
Borel space lower than the lower limit is considered invalid. As a measure of the relative 
reliability of various sum rules we consider the size of the regime in Borel space where both 
sides of the sum rules are valid. This quantity is used as a weight in the calculation of the 
I^DF as indicated in (|4.9| ). In addition, the size of continuum contributions throughout 
the Borel region can also serve as a measure of reliability, with small continuum model 
contributions being more reliable. 

The fiducial Borel region is chosen |p such that the highest-dimensional operator(s) 
(HDO) contribute no more than ~ 10% to the QCD side while the continuum contribution 
is less than ~ 50% of the total phenomenological side (i.e., the sum of the pole and the 
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continuum contribution). The former sets a criterion for the convergence of the OPE while 
the latter controls the continuum contribution. While the selection of 50% is obvious for 
pole dominance, the selection of 10% is a reasonably conservative criterion that has not 
failed in practice.^ See Ref. 0] for further specific examples. 

Figure ^ displays the valid Borel window for the p- meson sum rule of ( [5.1| ). The HDO 
contributions limited to 10% of the OPE and the continuum model contributions limited 
to 50% of the phenomenology are illustrated. The continuum contributions are the order 
of 15% of the phenomenological contributions at the lower end of the Borel regime. Here, 
the pole of interest truly dominates the QCD-SR. Figure ^ displays the corresponding fit 
obtained by minimizing in a three parameter search including /p, Mp and Wp. Figure 
§ displays the distribution of p-meson masses obtained from the 1000 sum rule fits. The 
distribution corresponds to 

= 0.76^|j;}^ GeV. (5.3) 

The uncertainty is somewhat larger than the commonly assumed 10%. Figure ^ displays the 
distribution of the square of fp which corresponds to 

/p = 0.144 ±0.021 GeV, (5.4) 

which agrees with the experimental value of fp = Mp/gp = 0.153 GeV. The continuum 
threshold is Wp = 1.23 ± 0.17 GeV which compares favorably with 1.31 ± 0.04 GeV (the p' 
mass less the half width), given the approximate nature of the continuum model. 



A. Derivative Sum Rules and the Ratio Method 



In addition to the fundamental sum rule of ( ^.1| ) one also finds the use of sum rules 
obtained by differentiating (|5.1| ) with respect to 1/M^. For the n'th derivative 



{MlY e-^l^'^ = n\ Co {M'r+' E„ 



Cn+2 
Cn+1 + TT^ + 



M2 



(m 



n 



1)! {M' 



2\m—n—l 



+ 



(5.5) 



While these sum rules are obviously unnecessary for the extraction of p-meson properties, 
the hope is that these additional sum rules will provide new useful information that might 
be used to further constrain the phenomenological properties. 

It is possible to isolate the p-meson mass as a function of the Borel mass, by taking a 
ratio of the first derivative sum rule with the fundamental sum rule. This approach has 
some aesthetic appeal, and has widely become the method of choice for analyzing QCD- 
SRs. The continuum threshold is selected to make the ratio of the two sum rules as fiat as 



^Reasonable alternatives to the 10% and 50% criteria are automatically explored in the Monte- 
Carlo error analysis, as the condensate values and the continuum threshold change in each sample. 
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FIG. 2. The valid Borel window for the p-meson sum rule of ( |5.1| ). The relative HDO contri- 
butions limited to 10% of the OPE and the continuum model contributions limited to 50% of the 
phenomenology are illustrated. 
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FIG. 3. A three parameter fit including /p, Mp and Wp of the /o-meson sum rule of ( p.l| ) for 
the central values of the condensate distributions. The error bars illustrate the OPE uncertainty 
estimated at each of the 51 points considered in the ^ estimate. Both the QCD— continuum 



(dashed) and pole (solid) sides of (l5.l|) are plotted 
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FIG. 4. Histogram of /9-meson masses obtained from 1000 fits of (|5.l| ) for 1000 QCD parameter 
sets. The distribution corresponds to Mp = 0.76 tgie GeV. 
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FIG. 5. Histogram of the squared p-meson decay constant obtained from 1000 fits of (|5.l| ) for 
1000 QCD parameter sets. The distribution corresponds to = 0.021 zt 0.006 GeV^. 
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possible as a function of Borel mass. Finally, the p-meson mass is selected from the point at 
which the ratio is most flat or stable. Unfortunately, this method has some very unpleasant 
shortcomings. 

The ratio method does not check the validity of each individual sum rule. It is possible 
to have individual sum rules that are not valid while their ratio is fiat as function of Borel 
mass. In addition, the ratio method cannot account for the fact that sum rules do not work 
equally well. The Borel region where a sum rule is valid can vary from one sum rule to 
another. Moreover, large uncertainties in one sum rule can spoil the predictions of another 
sum rule. Finally, the continuum contributions to the sum rules are not monitored in the 
ratio method. If the continuum contribution is dominant in a sum rule, one should not 
expect to get reliable information about the lowest resonance. 

The fact that sum rules do not work equally well deserves further discussion. In Ref. 
a detailed comparison of fundamental and derivative sum rules is carried out for the 
p meson. These arguments are not broadly known and are generally applicable to any 
derivative Borel-improved QCD-SR. Hence the important observations are highlighted here. 

If one could carry out the OPE to arbitrary accuracy and use a spectral density indepen- 
dent of models for excited state contributions, the predictions based on and those based 
on the derivative sum rules of ( |5.5| ) should be the same. In practical calculations, however, 
one has to truncate the OPE and use a simple phenomenological ansatz for the spectral 
density. Thus it is unrealistic to expect the sum rules to work equally well. Compare the 
n'th derivative sum rule of flS.SI) with the direct sum rule of (15.11). 



1. The perturbative contribution in the derivative sum rule has an extra factor n\ relative 
to the corresponding term in (|5.1|), implying that the perturbative contribution is 



more important in the derivative sum rules than in the fundamental sum rule, and 
becomes increasingly important as n increases. Since the perturbative term mainly 
contributes to the continuum of the spectral density, maintaining dominance of the 
lowest resonance pole in the sum rule will become increasingly difficult as n increases. 



2. In ( |5.1j ), the term proportional to Cm is suppressed by a factor of l/(m — 1)!, while it 
is only suppressed by l/(m — n — 1)! in the derivative sum rule (m > n). This implies 
that the convergence of the OPE is much slower in the derivative sum rule than in the 
fundamental sum rule. Consequently, the poorly known high-order power corrections 
are more important in the derivative sum rule than in Eq. ( |5.1| ), and become more and 
more important as n increases. If one would like to restrict the size of the last term of 
the OPE to maintain some promise of OPE convergence, the size of the Borel region 
in which the sum rules are believed to be valid is restricted. 

3. The power corrections proportional to Ci, C2, • • • , do not contribute to the n^^ deriva- 
tive sum rule but do contribute to ( |5.1| ). If one truncates the OPE, part or all of the 
nonperturbative information will be lost in the derivative sum rules. It is also worth 
noting that the leading power corrections are the most desirable terms to have. They 
do not give rise to a term in the continuum model and they are not the last term in the 
OPE, whose relative contribution should be restricted to maintain OPE convergence. 



In practice, the predictions based on the fundamental sum rule of (|5.1j ) are more reliable 
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than those from the derivative sum rules, which become less and less reliable as n increases. 
This can also be demonstrated by analyzing the sum rules numerically. 

Figure ^ displays the valid Borel windows for a simultaneous fit of the p-meson sum rules 
of i f) .11) and the first derivative [n = 1) of ( p.5|) . The reliability of the fundamental sum rule 
over the derivative sum rule is obvious. The valid Borel regime for the derivative sum rule 
is nearly nonexistent. 

All three points of the discussion surrounding derivative sum rules are displayed in this 
figure. Item 1 indicates the continuum model contributions will be much larger in the 
derivative sum rule for a given Borel mass. Indeed the continuum model contribution never 
drops below 45% in the valid regime. The relative contribution of the last term in the 
truncated OPE is larger in the derivative sum rule for a given Borel mass, in agreement 
with the expectations of items 2 and 3. 

Figure |^ displays the corresponding fit obtained by minimizing the weighted x^/^df of 
( [4.9| ) in a three parameter search including fp, Mp and Wp. This time the logarithms of the 
left- and right-hand sides of equations and ( |5.5| ) are plotted as a function of the inverse 



squared Borel mass. This method has the advantage of revealing the possible breakdown 
of the OPE signified by deviations of the dashed curves from linearity. Uncertainties are 
displayed only at the ends of the Borel region for clarity. 

While the central values of the best fit parameters Mp = 0.77 ± 0.07 GeV, = 0.022 ± 
0.002 GeV^, and Wp = 1.31 ± 0.05 GeV, remain essentially unchanged, the reduction of the 
uncertainties comes from rather unreliable information as discussed above. These arguments 
are further supported by the failure of the two sides of ( ^.5|) forn = 1 to agree as illustrated 
in figure |^. To further probe the reasons for this failure, the analysis of these sum rules is 
repeated with a spectral representation that accounts for the finite width of the p meson. 

To account for the finite width of the p-meson, the 6 function in ( [^.12D is replaced by a 
normalized Breit-Wigner form 

1 MpTp 

vr (s-M2 + r2/4)2 + M2r2 ' ^ ■ ' 

and the integral over the p-meson ground state 

is handled numerically Using the experimental width of 151.5 MeV, a three parameter fit 
of the fundamental and derivative sum rules with the Breit-Wigner form leads to essentially 
the same results reported in figure |^. The mass, coupling, and continuum threshold all 
agree within uncertainties, and have similar uncertainties. Hence, these spectral properties 
may be reliably estimated without precise information on the p width. However, p — uj 
mixing provides an interesting example where finite mesonic widths can play a central role 
in arriving at the correct physics [[71,^. 

The origin of the discrepancy for the derivative sum rule is more likely to lie in short- 
comings of continuum model. For this sum rule the continuum model contributions never 
drop below 45%, which may be too large for the simple model typically used in QCD-SR 
calculations. Alternatively, poorly known higher order terms, not present in the truncated 
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FIG. 6. The valid Borel windows for the fundamental p-meson sum rule of (dashed) and 
the first derivative sum rule of (|5.5| ) (dot-dashed). Both the relative HDO contributions limited to 
10% and the continuum model contributions limited to 50% are illustrated. 
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FIG. 7. A three parameter fit including /p, Mp and of the fundamental p-meson sum rule of 
(|5.1| ) (dashed) and the first derivative sum rule of ( |5.5D (dot-dashed) for the central values of the 
condensate distributions. Solid lines indicate the ground state contributions to these sum rules. In 
this and the following plots displaying sum rule fits, uncertainties are displayed only at the ends 
of the Borel region for clarity. 
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OPE may become important in the derivative sum rule. In the derivative sum rule, the next 
term in the series is suppressed by merely 2! . 

The importance of carefully monitoring the validity of QCD-SRs is nicely displayed 
in the recent debate over the behavior of vector-meson masses in nuclear matter. There 
an analysis involving invalid sum rules leads to the conclusion that vector-meson masses 
increase in finite density in contradiction to a previous analysis ||52| also based on 
the ratio method. The debate is resolved in Ref. [^] where valid sum rules indicate that 
vector-meson masses decrease for increasing nuclear matter density. 

Having firmly established the analysis procedure in the well understood p-meson channel, 
we are ready to apply the same approach to the nucleon sum rules. We will find that many 
results contradict the conventional wisdom. 



VI. NUCLEON SUM RULES: SPIN-1/2 INTERPOLATORS 



Here we focus on the Borel improved QCD Sum Rules for the generalized interpolator 
of ( p^.3| ). The sum rules proportional to the correlator Ils{M) where continuum model 
contributions are subtracted from the phenomenological side and placed on the OPE side 
are 



and 



5 + 2^+5^^6^-4/9 

64 



(1) 



256 

: a 

24 

— e , 



2M4 ^ M2 



^4/9 _ 13-2/3-11/3^ mja' 
96 M2 



27 



7-2(3- 5/3^ 
16 

16 



1 _ (.~n,yAI^ 



Wn 



M2 



19 + 10(3 - 29/32 



27 32 



a b 



(6.1a) 



(6.1b) 



where the Dirac-7 structures are indicated on the left, and = (2tt)'^Xc). The factors of 
L = [log(M/AQCD) / log (p/Aqcd)] account for the anomalous scaling of the operators and 
interpolators in the leading logarithmic approximation |]53|. 



A. Traditional QCD-SR Analysis 

Let us begin with the traditional or standard nucleon sum rule analysis where correc- 
tions are ignored. The more common interpolator of ( p.5| ) is utilized by setting /? = — 1 in 



27 



( |6.1|) . In addition, the traditional condensate parameter set with b = 0.47 ± 0.20 GeV^ is 
selected in order to evaluate the validity of previously published results. For a preliminary 
analysis of the validities of the sum rules, the nucleon mass is fixed to 1 GeV. Figure ^ illus- 
trates the valid Borel regimes for these two sum rules found by implementing the procedure 
discussed for the p-meson. 

A comparison of figure § with the analogous figure ^ for the p-meson suggests the nucleon 
sum rules will not be as reliable. In the nucleon sum rules, the continuum model contributes 
20% or more to the phenomenological side of the sum rules. In particular, the continuum 



model contribution in the sum rule of ( 6.1a ) exceeds that for the derivative p-meson sum 
rule which has already been shown to be incompatible with the fundamental sum rule. It 
should not be surprising to find a similar disaster for the nucleon case. 

As in the p-meson case, there is an analogous ratio method for obtaining the nucleon 
mass that has also become popular. Here the nucleon mass is isolated as a function of the 
Borel mass by taking the ratio of equations ( |6.1b|) / ( 6.1a|) . The mass is usually extracted from 



a regime in which the ratio is flat. Again OPE convergence and ground-state dominance 
criteria are not checked in this approach, and too often the mass is extracted from a regime 
in which the continuum model dominates. 

Figure || illustrates the three parameter fit of the residue Xq, and the two continuum 
model thresholds of wi and W2- The small Borel regime and the relatively huge OPE 
uncertainties of the sum rule of ( |6.1a| ) indicate that this sum rule plays a negligible role in 



determining the spectral properties of the nucleon. This is a perfect example of a case where 
large uncertainties in the sum rule of ( |S.la| ) would spoil the predictive ability of ( |6.1b| ) if 



the ratio method was used. The origin of the large uncertainties lies in the poorly known 
value of the four-quark condensate. The resultant fit parameters for this sum rule are also 
unphysical as the continuum threshold lies at zero below the nucleon pole. In other words, 
the present form of the sum rule at (7-^) is unable to provide any information on the ground 
state that is consistent with the more reliable sum rule of ( |6.1b| ) at the structure 1. 



Unfortunately, the community has incorrectly identified (|6.1a| ) as the more reliable sum 



rule. This belief is based on the simple fact that the leading term of the OPE (the identity 
operator) has a well known vacuum expectation value. However, the corrections for the 
identity operator are uncomfortably large. The first order corrections are independent 
of P and are 

""■''+7.), (6.2) 



TT V12 

where je = ^(1) — 0.58. This correction is ~ 50% and suggests perturbative corrections 
are completely out of control for this sum rule. Since the magnitude of the identity operator 
contribution is closely tied to the size of continuum model contributions, it is not surprising 
to find the valid Borel regime for this sum rule disappears when corrections are added. 
Another factor in the traditional preference of (|6.1a|) over (|6.1b|) is that the second term 



of ( |6.1b| ) vanishes for loffe's choice of nucleon interpolator [jS = —1). This is really nothing 
more than a good reason to exclude the interpolator xsr of (|2.5| ) as an optimal interpolator. 

When one considers the monstrous leading-order corrections, the possibility of a sig- 
nificant dimension-two power correction arising from a summation of the perturbative series 



25,21], and the large uncertainties associated with the value of the four-quark condensate, 
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FIG. 8. The valid Borel windows for the nucleon sum rules of ( 6.1a| ) (dashed) and ( |6.1b ) 
(dot-dashed). Both the relative HDO contributions limited to 10% of the OPE and the continuum 
model contributions limited to 50% of the phenomenology are illustrated. 
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FIG. 9. Three parameter fit of (|6.1a| ) (dashed) and (|6.1b|) (dot-dashed) to the ground state 
contribution (solid) for the residue Xq, and the two continuum model thresholds of wi and W2- 
The small Borel regime and the relatively huge OPE uncertainties of the sum rule of ( 6.1a| ) indicate 
that this sum rule plays a negligible role in determining the spectral properties of the nucleon. 
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one cannot help but reach the conclusion that the QCD side of the sum rule of ( |6.1a|) is 
really unknown. 



B. Nucleon Sum Rule at 7-p 



Since the sum rule of ( |6.1a ) has traditionally been the nucleon sum rule of choice, it is 



interesting to consider the predictions of this sum rule alone. Applying the same techniques 
as used for the p-meson to ( |6.1a| ), a three parameter fit of A^, Wi and yields 

= 0.15 ±0.05 GeV^ Wi = 0.74 ± 0.19 GeV , Mjv = 0.44 ± 0.05 GeV . (6.3) 

While setting kat = 1 slightly increases the nucleon mass to 0.47 GeV, adding leading order 
Us corrections has a negligible effect on the fit parameters. Consideration of the more recent 
condensate parameter set discussed in Sec. fT^ leads to similar results 

A^ = 0.14 ±0.10 GeV^ Wi = 0.71 ± 0.5 GeV , = 0.43 ± 0.09 GeV . (6.4) 

This may be the first QCD sum rule to predict spectral properties that significantly contra- 
dict experiment .[] As such, this may be the most interesting of all QCD sum rules, and we 
shall return to this sum rule after examining the selection of the most optimal interpolator 
for nucleon sum rule analyses. 

The predictions of (|6.1a| ) summarized in (|6.3|) are extremely different from that found in 



the literature. Hence it is important to understand how these new results come about. A 
critical examination of the nucleon ground state results reported in Ref. by the present 
author sheds considerable light on the importance of the rigorous QCD-SR analysis presented 
here. 

Ref. was the first analysis to give careful regard to the issues of OPE convergence 
and ground-state dominance of the phenomenological side. Moreover, all phenomenological 
parameters were included as fit parameters to be determined by the sum rules. However, 
there are three important differences between the analysis presented here and that of Ref. 



The standard logarithmic measure |^ of sum rule agreement was used to fit the sum 
rules. This approach fails to apply an appropriate weight to sum rule discrepancies, 
unlike the measure of 



2. The algorithm used to fit the sum rules did not employ the cautious convergence 
criteria discussed in Section [IV B| . Implementation of the new criteria aids in avoiding 
insignificant local minima. 

3. Only the standard condensate parameter set was considered in determining the nucleon 
mass. Here 1000 sets are considered. 



^Note that this is unlike the case of the pseudoscalar pion current |55| where a check of OPE 
convergence and ground state dominance indicates the absence of a valid Borel regime. 
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The optimal parameter set found in Ref. corresponds to a local minimum which is sensitive 
to the measure of sum rule agreement, the initial parameter estimates and the condensate 
parameter set. 

Other authors have managed to keep the nucleon mass closer to the experimental result 
by extracting the result from a regime in Borel space that provides the correct answer or 
fixing the continuum threshold to values suggested by the particle data tables. As we shall 
see in Section |V111| , fixing the continuum threshold to phenomenologically motivated values 
largely determines the mass. Hence you can get anything you want provided you are clever 
enough to fix the threshold to the right value. And since the excited states are modeled 
there is tremendous freedom in the threshold parameter space. Of course, such input may 
not be in agreement with the content of the OPE. In any event, such results must not be 
regarded as predictions of QCD sum rules. 



C. Optimal Mixing of Spin-1/2 Interpolators 



Having established the nucleon sum rule of ( |6.1bD at the structure 1 to be the poten 



tially more reliable sum rule, it is important to determine the optimal mixing of spin-1/2 
interpolators for this sum rule. 

It has long been established that there are two independent interpolating fields with no 
derivatives having the quantum numbers of spin 1/2 and isospin 1/2. Both are expected to 
excite the ground state nucleon from the vacuum. As such, there is a history of argument 
over the optimal nucleon interpolating field to be used in analyses [pi |T0tp!2| - [T5|j56| -^ and 
this issue is recognized to be of paramount importance [^ . 

While some advocate an interpolating field for which the leading terms of the OPE are 



stationary with respect to the interpolating field mixing parameter []T0|,^, others argue that 
a balance between OPE convergence and QCD continuum contributions must be maintained 
P JT5|j56|j57|| . Ideally, one would like to simply calculate with alternate interpolating fields 



and confirm that the nucleon properties remain unchanged. However, previous implemen- 
tations of the QCD-SR approach have hindered any firm conclusions. The rigorous analysis 
introduced here will shed considerable light on these issues. 

Fortunately, the lattice approach to QCD is not plagued with the same limitations and 
the debate over the optimal interpolating field was resolved in Ref. [Q. There is was 
found that, to a good approximation, X2 of ( |2.4b| ) excites pure QCD continuum. Since 
X2 has negligible overlap with the ground state nucleon, it is tempting to simply conclude 
that the optimal interpolating field is Xi- While this is certainly the case for lattice QCD 
investigations, it is not the case for QCD-SR analyses. By including X2 components in an 
interpolating field, one can reduce the continuum contributions excited by xi and allow a 
broader Borel analysis window in which better contact with the ground state of interest is 
obtained. 

One of the most difficult things to monitor in the QCD-SR approach is whether the 
OPE is sufficiently convergent for a particular value of Borel mass. The lattice results of 
Ref. [0 indicate the X2X2 correlator has the fastest converging OPE, as its overlap with 



the nucleon ground state is negligible. Similarly, the combination XiXi produces an OPE 
with the slowest convergence, as this correlation function is dominated by the ground state 
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FIG. 10. Continuum model (solid curve) and HDO (dashed curve) contributions to the Borel 
improved QCD-SR of ( 6.1b| ) plotted as a function of the interpolating field mixing parameter p. 



nucleoli for small Borel masses. 

Hence, errors made in truncating the OPE are dominated by errors in the XiXi com- 
ponent of the general correlator. The relative error in the OPE truncation can be reduced 
by adding X2 components to the correlator. However, the X2 components in the OPE are 
simply subtracted off again by the continuum model terms. Hence the relevant error is the 
absolute error. For \f3\ < 1, this error is dominated by XiXi components of the correlator. 
As a result, OPE truncation errors are approximately independent of (3. This crucial point 
has been neglected in previous arguments regarding the optimal nucleon interpolating field. 

Since the errors made in truncating the OPE are dominated by errors in the XiXi com- 
ponent of the general correlator the criteria of a 10% limit on the HDO contribution to the 
OPE should be applied to sum rules when j3 = 0. Once this lower limit is established, it 
may be used for all < 1. We will refer to this approach for determining the valid Borel 
regime as the "refined" method in the following. 

The optimal nucleon interpolator must involve Xi, as this interpolator is required to 
maintain overlap with the ground state. The task is to determine the optimal mixing of 
X2- Since X2 has negligible overlap with the nucleon, the ground state contribution to a 
sum rule is also independent of /3. Hence, the size of the continuum model contribution is 
the predominant factor in determining the optimal interpolator. Figure |1^ illustrates the 
contributions of the continuum model terms in (|6.1bD for M = 0.94 GeV and W2 = 1.4 
GeV. The following discussion is not dependent on the precise values of these parameters. 
The first point to be made is that contributions from the continuum model are largest for 
f3 ~ —0.2. This selection of mixing is the worst possible choice for extracting information 
on the ground state nucleon. 
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Figure ^ also indicates it is possible to have vanishing continuum model contributions at 
1.5or/5 = l. However, we are relying on the continuum model to account for strength 
in the correlator that does not have its origin in the ground state. Without a continuum 
model, one would need to include additional poles on the right-hand-side of ( |6.1b|) to account 
for positive and negative parity excitation strength. For (3 < —1.5 or /5 > 1.0 the correlator 
is negative indicating the sum rule is saturated by a negative parity state. 

Thus the optimal interpolator is /5 ~ —1.2 or /? ~ 0.8. To discriminate between these two 
regimes, we turn to the higher-dimension-operator contributions which do not contribute to 
the continuum model. It is these terms that provide crucial information on whether the 
strength in the correlator lies in the ground state or the excited states. If these terms are 
absent, the optimal fit of the correlator is obtained when Ao — and W2 0. In this 
case the continuum model and the OPE become equivalent via the Laplace transform of 
( p.26| ) and the fit is perfect. Hence the higher-dimension-operator contributions should be 
large in magnitude. A change in sign from the leading terms of the OPE will also assist 
in distinguishing ground state strength from excited state strength as the change in the 
curvature of the correlator will be more prominent. 

The last term of (|6.1b|) is independent of the continuum model, and its value is plotted 



as a function of beta in Figure 0. The contributions are larger for (3 ~ —1.2 than for 
f3 ~ -1-0.8. In addition, the sign of the contribution is opposite that of the continuum model 
contributions. Hence the preferred regime is (3 ~ —1.2. These arguments are supported by 
the numerical analysis described below. 

In summary, the lattice results indicate OPE convergence, and ground state pole contri- 
butions are approximately independent of /3. Considerations of the size of continuum model 
contributions and the sign and magnitude of operators independent of the continuum model 
leads to the preferred value of 

/5 = -1.2 ±0.1. (6.5) 

A more precise determination of (3 will depend on the details of limits for continuum model 
and HDO contributions, condensate values, and other parameters of the sum rules. 



At (3 = —0.2, where the leading term of the OPE is stationary with respect to (3 p0| , |58 
the continuum contributions are maximal. The positive value and small magnitude of the 
HDO contribution indicates the stability of the leading terms of the OPE will not be realized 
as stability in the ground state mass, coupling, nor in the continuum threshold. 

Another point of view argues that the optimal interpolating field is one that eliminates 
direct instanton contributions from the correlator However, such arguments are based 



on the assumption that direct instantons necessarily contribute to the nucleon correlator in 
a significant manner. Section |^ will demonstrate that direct instantons are not necessary 
to maintain sum rule self consistency for a wide range of interpolating fields. The f3 depen- 
dence of the direct instanton contributions to ( |6.1b| ) is proportional to 13/5^ -|- lOP -|- 13 7^ 0. 
Minimization of instanton contributions would place f3 = —5/13, very near (3 = —0.6 where 
the HDO term of ( |6.1b ) vanishes and where continuum contributions are near their peak. 



Isolation of ground state contributions is nearly impossible in this regime. Therefore mini- 
mization of direct instanton contributions does not lead to an optimal interpolating field. 

While it is important to establish the optimal mixing of interpolating fields for QCD-SR 
analyses, one should not overlook the fact that there is a range of values for /3 where the 
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sum rules are expected to work. Moreover, the ground state contribution to all these sum 
rules is equivalent to the 1% level [0. In other words, the right-hand side of ( |6.1b| ) for a 
single pole plus continuum model is independent of (3. After the first sum rule is written 
down, additional sum rules may be introduced with merely one new fit parameter (the 
continuum threshold) per sum rule. Since direct instanton contributions to the sum rules 
are not independent of /5 0], one has an excellent opportunity to see if direct instanton 
contributions really are necessary to maintain sum-rule consistency. This issue is addressed 
in Section |X| below. 



D. "loffe Formula" 



Having discovered that X2 has negligible overlap with the nucleon ground state, it might 
be interesting to return to the "loffe Formula" this time generalized for arbitrary j3 



M, 



N 



7-2/3-5/32 
5 + 2/3 + 5/32 



4 {2ny {qq) 
M2 



(6.6) 



For p = 0, Mn = 2.5 GeV when M = 1 GeV, or Mn = 1.4 GeV when M = Mn GeV. 
The explicit dependence on the choice of interpolating field and the Borel mass limits the 
significance of the "loffe formula" to a qualitative role at best. 



E. Nucleon Sum Rule at 1 for 0=-1.2 



Here we focus on the more reliable sum rule of ( |6.1b|) at the optimal interpolator mixing 



of /3 = —1.2. Fortunately, the corrections for the leading order operator of the OPE are 
much smaller than for the identity operator at the order of 10%. The /3-dependent coefficient 
of the quark condensate in ( |6.1b| ) becomes 



These corrections are small relative to uncertainties in the other terms of the correlator and 
can be absorbed by a small shift in the condensate itself. Since our QCD parameter set is 
extracted from analyses with neglected as corrections, we shall do the same here. 

Armed with knowledge of the optimal interpolator, and a refined understanding of OPE 
convergence tests, we now present the QCD-SR predictions for the nucleon spectral prop- 
erties from spin-1/2 interpolators. Figure |TT| displays the valid Borel regime for this sum 
rule. The HDO contributions are less than 10% of the OPE contributions for this fit when 
/3 = 0, as indicated by the solid curve. As described above, contributions from correlator 
components known to be more convergent increase the HDO operator contributions making 
resolution of the pole from the continuum more probable. 

Figure |T2| displays the fit for this sum rule. The near perfect linearity of the 
OPE— continuum side of Ils{M) supports the arguments presented above. Despite 40% 
contributions of the HDO term to this sum rule, the OPE— continuum side displays the 
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FIG. 11. The valid Borel regime for the sum rule of ( |6.1b[ ) at the optimal f3 = —1.2. Relative 
HDO and continuum model contributions are illustrated (dot-dash curves). The HDO contributions 
are less than 10% of the OPE contributions for this fit when /3 = 0, as indicated by the solid curve. 
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FIG. 12. Three parameter fit of ( |6.1aD for the residue Xq, the continuum model threshold W2, 
and the nucleon mass. The near perfect linearity of the OPE— continuum side (hidden dot-dash 
line) of 115 (M) supports the arguments surrounding the refined determination of the valid Borel 



regime . 
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anticipated behavior of reasonable OPE convergence. The optimal fit parameters obtained 
from consideration of 1000 QCD parameter sets are 



A 



o 



0.20 ±0.12 GeV^ 



W2 



1.53 ± 0.41 GeV, Mat = 1.17 ± 0.26 GeV . 



(6.8) 



The uncertainty of the nucleon mass is more than double the traditional expectation at 260 
MeV. Figures 0, |T^, and |TB| display histograms for the nucleon mass, threshold and residue 
respectively. These histograms will aid in the discussion of the following section. 

Table |I| compares the prediction for Xq with other estimates of Asr = 2 Xq- This residue 
corresponds to loffe's interpolator of and is expected to be independent of < 1 to 



the order of 1% |3| 



In the past the community has emphasized the importance of the independence of the 
nucleon mass from the Borel mass parameter. It is common to find a plot of the nucleon 
mass extracted from the ratio method discussed above as a function of the Borel mass. While 
independence of the spectral parameters from the Borel mass is certainly desirable, the ratio 
method is potentially problematic in practice. We have already seen that the reliabilities 
and validities of two sum rules are usually different and could even result in separate valid 
regimes. This feature cannot be revealed in the ratio of the two sum rules. In addition, one 
can always achieve flatness in the ratio method in the large Borel mass region, where both 
sides of the sum rules are dominated by the continuum model. However, one learns little 
about the lowest pole in this case. 

Contact with the plateau criterion may be made by rewriting the sum rule of interest 
solving for the spectral parameter in question. For example, one can express the nucleon 
mass in the sum rule of ( |6.1b|) as. 



Mm = M 



In 



)i2 



.1/2 



(6.9) 



where Ils{M) denotes the left-hand side of (|6.1b|) . Indeed a plot of versus the Borel 
parameter M is very fiat within the valid region. This information is already apparent 



TABLE I. Comparison of predictions for Asr = 2 Ao for various approaches to QCD. 



Approach 


Reference 


Asr 






(xlO-2 QgV3) 



QCD Sum Rule 
Lattice (mean-field improved) 
Lattice (conventional renormalization) 
Lattice (coordinate space) 
Instanton Liquid 
Baryon wave functions (x^ 
Quark Model 
Bethe-Salpeter amplitude 
Quark Model 
MIT Bag Model 



0) 



This work 
Leinweber []l3| 
Gavela et al. |61| 



Chu et al. |2| 
Schafer et al. 



Brodsky et al. |64| 
Thomas and McKellar 
Tomozawa 



Milosevic et al. §7\ 
Donoghue and Golowich 



2.2±0.7 
2.7±0.5 
2.4 

2.2±0.4 
3.2±0.1 
12 
8 

2.5 

2 

1.27 
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FIG. 15. Histogram for the residue of the nucleon pole obtained from fits of ( |6.1b ) at the 
optimal (3 = —1.2 for 1000 QCD parameter sets. 



in figure |T2] where the near perfect Unearity of the OPE— continuum side of the sum rule 
indicates independence of the fitted mass from the Borel parameter. For a similar analysis 
of the p-meson in finite density nuclear matter see . 

Returning to the sum rule at the structure 7 ■ p of ( |6.1a| ), one finds the valid Borel 
regime disappears when the refined OPE convergence criteria is applied. The dashed curves 
of figure |16| display the naive Borel regime for ( |6.1a| ) at /? = — 1 used in the fit of ( |6.4|) . 
The solid curve illustrates the HDO contributions to the slowest converging component of 
the correlator [f5 = 0). This curve fails to drop below 10% before the continuum model 
contributions exceed 50%. The HDO of the faster converging components reduces the HDO 
contribution of the slowest converging component to give the illusion of reasonable OPE 
convergence at (3 = —1. This is most likely the explanation for the failure of the sum rule 
of ( |6.1aD ; the absence of a valid Borel regime. 



F. Nucleon Sum Rule at 1 for 0= +0.8 

Having considered the most favorable interpolator for the nucleon, let us now examine 



the QCD-SR predictions for a less than favorable interpolator. At /5 = +0.8, figure ^ 
indicates the continuum model contributions can be reasonably small. Hence it should be 
possible to extract ground state information that is reliable. However, the HDO term is of 
the same sign and small relative to the leading order contributions. The difficulty here will 
be resolving the ground-state pole from the approximate continuum model. 

The lower limit of the valid Borel regime is determined by the OPE contributions which 
are independent of the phenomenological fit parameters and independent of /?. Hence this 
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FIG. 16. The naive Borel regime for ( 6.1ei| ) at /? = —1. The sohd curve ihustrates the HDO 
contributions to the slowest converging component of the correlator (/3 = 0). The failure of this 
curve to drop below 10% before the continuum model contributions exceed 50% indicates the 
absence of a valid Borel regime. 



limit is known to be 0.88 GeV. The inherent instability of this sum rule makes it difficult 
to determine the upper limit. Figure |T0| suggests that continuum contributions should be 
similar to that for j3 = —1.2 and therefore we take an upper limit of 1.17 GeV, as in Figure 
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Figures |17|, [1^, and ^ display histograms for the nucleon mass, threshold and residue 
obtained from ( 6.1b ) for j3 = +0.8. Figure T7 resolves two regimes for the nucleon mass and 
figure suggests the presence of at least three regimes. 

In section |V111] the nature of these three regimes are resolved via a correlation analysis. 
As mentioned earlier there are two ways to fit a sum rule. If the terms which do not 
contribute to the continuum model are insufficient to resolve the pole, a good fit may be 
obtained by letting the continuum threshold and the residue of the pole go to zero. One 
then has an approximate Laplace transform relating the two sides of the sum rule. These 
fits are represented by the left-most shoulder of the histograms. 

An intermediate regime arises when the sum rules are capable of resolving a threshold at 
which strength in the spectral density begins. However, a gap between ground and excited 
states is not resolved. The peak at 0.5 GeV in the continuum threshold is representative of 
these fits. For these fits the nucleon mass is scattered about typical values but with with a 
small residue. The smaller residue compensates for strength from the continuum model. 

Finally, the third regime represented by the broad right-most peak in the continuum 
threshold histogram corresponds to more standard fits where a gap between the ground and 
first excited state has been resolved. 
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FIG. 17. Histogram for the nucleon mass obtained from fits of ( |6.1b ) at (3 = +0.8 for 1000 QCD 
parameter sets. When the pole is resolved from the continuum, the histogram compares favorably 



with that of figure 13 




1.0 1.5 
(GeV) 



2.5 



FIG. 18. Histogram for the continuum threshold obtained from fits of (6.1b) at /3 = +0.8 for 
1000 QCD parameter sets. The three types of sum rule fits are clearly resolved. 
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FIG. 19. Histogram for the residue of the nucleon pole obtained from fits of ( |6.1b| ) at /3 = +0.8 
for 1000 QCD parameter sets. Additional counts at small allows for additional strength from 
the continuum model. 



It is encouraging to see the similarity in the nucleon mass distributions of figures ^ 



and |T3| when the sum rules resolve a pole plus continuum fit. This agreement comes about 
without the need for direct instanton contributions. 



VII. NUCLEON SUM RULES: SPIN-1/2 AND 3/2 INTERPOLATORS 

A. Sum Rule Derivation 



The spin-3/2 interpolating field for the nucleon is 



abc 



{u''^{x)Capxd\x)) a^'^^^u^x) - {u''^{x)Capxu\x)) aP^-f^'d'ix) . (7.1) 

Nucleon to vacuum matrix elements of X3/2 defined as 

(0|x^/2|l/2+) = {ap^ + Xy2Y)l5u{p), (7.2a) 
(0|x^/2|l/2-) = {ap^ + Xs/2r) u{p) . (7.2b) 

The appearance of 75 on the right-hand side matches the parity of Xz/2- "^^^ property 
7^ X3/2 — provides the relation 



4A 



a 



3/2 



±M. 



N 



(7.3) 
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In an analogous manner we work with the generahzed spin-1/2 interpolator 



Xo = 1^15 Xo , 



(7.4) 



which has the matrix elements of ( |7.2| ) with a = and A3/2 — ^ ^o- Overlap of the spin- 
3/2 interpolating field with the generalized spin-1/2 interpolator of (|2.3|) gives rise to four 
Dirac-7 structures 

(0|x^|l/2±)(l/2±|x3^/2|0) = AaA3/2 {-7V (7 ■ P ± M^) ± -^Yp'l ■ P - ^Yp"} , (7-5) 
from which three independent sum rules are obtained. 



(7^7, 7 -p) : KMa^L^'/^' 



9 — (3 mla? 
24 M2 



^2/9 



^0^3/2 e ^' 



(7.6a) 



2 




1 + 3/j 
96 



(7.6b) 



2 



3-/3 
16 



1 + 3/? g & ^8/27 

96 M2 



(7.6c) 



The corresponding Dirac-7 structure from which the sum rule is obtained is indicated on 
the left. The sum rule at the structure ■jfj.Pu is identical to that for ■jfj.'y^'y ■ P- 

The leading term of the sum rule at the structure 7;_t7jy 7 ■ p is the four-quark condensate. 
As such the OPE of this sum rule is poorly known. In addition, there are no terms from 
which a continuum model may be constructed. For an analysis and discussion of this sum 
rule under the assumption of factorization see 0. 

The sum rule at the structure 'y^Qu 7 • p is a very favorable sum rule. The dimensionality 
of the structure at which the sum rule is extracted allows the HDO terms of the sum rule 
to approach the factorial suppression regime of the Borel transform thus improving the 
truncated OPE convergence. Moreover, there are two terms which are not used in the 
continuum model. Hence there should be better resolution of the pole from the continuum. 

The mixed condensate is absent in ( [7.6b| ) for all j3. With only two pieces of information 
on the OPE side of the sum rule, this sum rule is of little use on it's own. Fortunately, 
both ( |7.6bD and ( |7.6c| ) are extracted from Dirac 7-structures in which even- and odd-parity 
excitations contribute with opposite signs. Hence it is reasonable to set the continuum 
model thresholds equal for these two sum rules and this provides greater stability in the fit 
parameters. 
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Since the terms of ( |7.6|) used in the continuum model are independent of the interpolator 
mixing parameter (3, the obvious selection for the optimal mixing is /? = 0. Of course, if 
the interpolator X2 of ( |2.4b| ) doesn't contribute to the continuum model, nor to the ground 
state, it will be important to demonstrate that different values of j3 lead to similar results. 
The expectation is that the continuum threshold will be dependent on (3 while the mass 
and residue remain unchanged for \(3\ < 1. This dependence is examined in Section |X 
examining the necessity of direct instanton contributions. 

loffe's selection of j3 = —1 as the optimal spin-1/2 interpolator is reasonably close to the 
determination of /? = —1.2 in Section [Vl (J| for spin-1/2 to spin-1/2 interpolators. However, 
the optimal mixing is specific to each individual sum rule. Here we have found an optimal 
value quite different from loffe's preference and, as we shall see, loffe's selection of /? ~ — 1 
fails to provide a sum rule with a valid regime. 



B. Sum Rule Analysis 



We begin by considering a simultaneous analysis of (|6.1b| ) at P = —1.2, and (|7.6b|) and 
( [7.6c| ) at /3 = 0. A valid Borel regime could not be found for ( |6.1b[ ) when analyzed in 
conjunction with the two additional sum rules of (|7.6b| ) and (|7.6c|) . Hence we will analyze 
( [7.6b| ) and ( [7.6cD alone. Figure |2y displays the valid Borel regimes for these two sum rules. 



The corresponding three parameter fit of the residue Xo^3/2, threshold, 1^3, and nucleon 



mass is indicated in figure The sum rule of (|6.1b|) is also displayed for the Borel regime 
considered in Sec. |V1E[ . The linearity of the OPE— continuum sides of the subtracted sum 
rules over a relatively broad Borel regime indicates the OPE is under reasonable control. 
However, since the continuum model dominates the phenomenology of ( |6.1b| ), we do not 
consider it further. The resultant fit parameters obtained from the consideration of 1000 
QCD parameters sets are 



M, 



N 



0.96 ±0.08 GeV. 



(7.7a) 



Ao A3/2 = 0.41 ± 0.14 GeV^ , = 1.3 ± 0.2 GeV . 



(7.7b) 



The uncertainty of 100 MeV for the nucleon mass is a significant improvement over the 
traditional sum rule of ( 6.1b ) alone as illustrated in figure The value of the threshold 
lies roughly at the value suggested by the particle data tables of approximately 1.27 GeV; 
the Roper mass of 1.440 GeV less the half width of 0.175 GeV. The value for A3/2 is a 
prediction that awaits vindication from alternative approaches. 

The presence of a valid Borel regime common to both sum rules of (|7.6b|) and ( [7.6c|) 
provides an opportunity to apply the Monte-Carlo uncertainty analysis to the ratio method. 
The ratio of equations ( [7.6b|) / (|7.6c|) provides a ground state contribution equal to and 
independent of the residues \oX3/2- The optimal fit parameters are 



M, 



N 



1.05 ±0.08 GeV, W3 = 1.7 ± 0.5 GeV , 



(7. 



and agree with the previous results. The difference in the central values is due to the 
different Borel regimes used in the fit. Despite taking a ratio of sum rules, the uncertainty 
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FIG. 20. The valid Borel windows for the nucleon sum rules of ( [7.6bD (dashed) and ( 7.6c ) 
(short dash-dot) at the optimal /? = 0. The relative HDO contributions limited to 10% and the 
continuum model contributions limited to 50% are illustrated. 
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FIG. 21. A three parameter fit including A3/2, Mn and ws of the nucleon sum rules of ( 7.6b ) 
(dashed) and ( 7.6c| ) (short dash-dot) at /? = 0. The sum rule of ( |6.1b| ) at /? = — 1.2 (long dash-dot) 
is also illustrated to demonstrate the consistency of all three sum rules. The solid lines illustrate 
the ground state contributions. 
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FIG. 22. Histogram for the nucleon mass obtained from fits of ( 7.6b| ) and ( |7.6c ) at the optimal 
/3 = for 1000 QCD parameter sets. 



in the nucleon mass remains nearly unchanged. The fit obtained from the full covariance 
matrix of provides similar results and yields x^/^df = 1-34. This acceptable value 
suggests the distributions selected for the QCD parameters are in reasonable accord with the 
approximations inherent in the QCD-SR approach. Future studies might use the x'^/^df 
to place restrictions on QCD parameter sets selected in the Monte-Carlo approach. 

Table H summarizes an investigation of x^/^df and uncertainties obtained from the 
analysis of ( [7.6 b| ) and ( [7.6cD at /5 = for various relative errors of condensate input param- 
eters. Results are based on samples of 200 condensate parameter sets. The ratio results are 
obtained from the ratio of these sum rules analyzed in the Borel regime common to both 
sum rules. A comparison of the uncertainties for the nucleon mass obtained from the ratio 
method and the standard analysis advocated here indicates there is very little if anything 
to be gained in the ratio approach to sum rule analysis. 

Table || indicates a x'^/^bf — 1 is obtained when the QCD input parameters are specified 
with 12% uncertainties. Hence, the best one can determine phenomenological parameters 
with the current implementation of QCD-SRs is 10% accuracy for the nucleon mass and 25% 
accuracy for the residue. This large uncertainty for the residue does not bode well for the 
analysis of three-point functions where the quantity of interest appears multiplied by the 
reside. The two-point function is used to normalize the three-point function, and it appears 
that this normalization is not well determined. Reduction of this uncertainty to the 10% 
level requires 5% uncertainties on the QCD input parameters. However, the x^/A^df = 18.7 
indicates this level of accuracy is beyond the current state of the art. Of course, this is not 
surprising since approximations at the level of 10% have been made in deriving the sum 
rules. 
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TABLE II. An investigation of x^/^df and uncertainties obtained from the analysis of ( [7.6b ) 
and ( [7.6c) ) at /3 = for various relative errors of condensate input parameters. Results are based 
on samples of 200 condensate parameter sets. Ratio results are obtained from the ratio of these 
sum rules analyzed in the Borel regime common to both sum rules. The last row summarizes the 



results for the condensate parameter set estimated in Sec. Ill 







Sum Rule Ratio 






Standard Analysis 




Condensate 




Mn 


W3 


Mn 




-^0^3/2 


Uncertainty 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV^) 


20% 


0.11 


1.04(12) 


1.78(91) 


0.96(11) 


1.29(26) 


0.40(15) 


15% 


0.26 


1.04(10) 


1.73(75) 


0.96(10) 


1.29(22) 


0.40(12) 


13% 


0.69 


1.04(9) 


1.72(68) 


0.95(9) 


1.28(19) 


0.40(10) 


12% 


1.13 


1.04(8) 


1.70(65) 


0.95(8) 


1.27(17) 


0.39(9) 


10% 


2.22 


1.04(6) 


1.67(47) 


0.95(7) 


1.26(14) 


0.39(8) 


5% 


18.74 


1.02(3) 


1.47(6) 


0.95(3) 


1.26(7) 


0.38(4) 


Estimated 


1.34 


1.05(8) 


1.72(52) 


0.96(8) 


1.28(20) 


0.41(14) 



VIII. CORRELATIONS 
A. Measures of Association 

To quantitatively examine the correlations among the QCD and phenomenological pa- 
rameters, we utilize a contingency table analysis of two distributions to determine the 
significance, and strength of correlations. For the correlations found to be significant, the 
linear correlation coefficient is calculated. 

Determinations of the significance, probability and strength of correlations can be based 
on measures of association. Null-hypothesis probabilities, Cramer's V, and the Contin- 
gency Coefficient C are a few of the better known measures of association which may be 
obtained from the contingency table analysis described below. 

A contingency table analysis of association between two variables proceeds by binning 
the two-dimensional distribution into an J x J grid [|^. If Nij is the number of points in 
bin i,j, then the row and column totals are 

N,. = j2Nij, N., = j2N^j, (8.1) 

j i 

and the total number of points in the distribution is 

^ = E% = E^^- = E^.- (8-2) 

The null hypothesis is the assumption of no association between the two variables. By 
measuring the fit of this assumption to the data we can learn about the significance of the 
association. Under the null hypothesis, the fraction of points in bin i,j, of column j is 
independent of j. That is. 
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where riij is the number of points in bin i,j predicted by the null hypothesis. The is 

= E ~ ""''^^ , (8.4) 
i,j "-u 

and the number of degrees of freedom are 

u = IJ-I-J+l, (8.5) 

where the row and column totals used in the hypothesis are subtracted and the addition of 
one accounts for the equality of the sum of row and column totals. 

Knowledge of x^/^ allows an estimate of the probability for the null hypothesis to be 
true. A small probability indicates significant association between the data. We calculate 
the standard null-hypothesis significance parameter 



NH 



T{iy/2) 7xV2 



oo 



-'f'/^-'dt, (8.6) 



where F is the gamma function. 

Cramer's ^ is a normalized measure of the strength of association independent of the 
size of the contingency table 

_ / \ ^'^ 

^ ~ \N Ymn{I - I, J - I)) ■ ^^'^"^ 

This measure lies in the range 

0<1^<1. (8.8) 

V = 1 indicates there is a unique square for each i,j holding all points in row i and column 
j. V = indicates there is no association. 

The Contingency Coefficient C is also a measure of the strength of the correlation. It is 
defined as 



and is dependent on the size of the contingency table. Since all the contingency tables used 
in this analysis are the same size, this poses no problem. Its range is 

0<C<1, (8.10) 

with larger values indicating stronger association. 

We also consider the significance of association based on the symmetrical uncertainty 
coefficient U determined via an examination of the entropy of the distribution P9|. In 
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essence this measure describes how much entropy is lost in one parameter given knowledge 
of the other. This measure is bounded by 



0<U<1, (8.11) 

and large values of U indicate significant loss of entropy implying significant dependency 
among the parameters. 

Once a distribution has been found to be significant, the linear correlation coefficient is 
calculated 



T,t{xi - x){yi - y) 



(8.12) 



The linear correlation coefficient also measures the strength of the association and lies in 
the range 

-l<r<l. (8.13) 

When r = +1 the points lie on a perfect line with positive slope, while r = —1 indicates 
the points lie on a perfect line with negative slope, r = indicates there is no correlation. 
While the linear correlation coefficient is a poor measure of the significance of correlation, 
it has the advantage of identifying the sign of significant correlations. 

Hence the procedure is as follows. The goodness of fit parameter is utilized to iden- 
tify associations that are significant. This parameter has the advantage of presenting a 
clean determination of which associations are significant. Cramer's V , the contingency co- 
efficient C, and the symmetrical uncertainty coefficient are used to determine which QCD 
parameters play the dominant role in determining hadronic spectral properties. Finally, the 
proportionality of the association is measured by the linear correlation coefficient. Since the 
uncertainties of the QCD parameter distributions can affect the strengths of the associations, 
we utilize 1000 QCD parameter sets in which all QCD parameters have an equal relative 
uncertainty of 15%. However, the illustrations of this section show the actual correlations 
for the QCD parameters described in Section |ITJ 



B. p-meson Associations 



Table |ITI| identifies the QCD parameters having significant correlations with the spectral 
properties of the p-meson. Here the quark condensate plays a dominant role in determining 
all the spectral properties. This is simply a reflection of the fact that it appears squared in 



the last term of (|5.1|) . The correlations for show a linearized version of essentially the 



same correlation. The a<j correlations are similar to that for Kp and reflect the fact that the 
Us corrections to the identity operator are small relative to one. An interesting aspect is the 
absence of significant gluon condensate correlations for and Wp. 

The correlations between the p-meson mass and two of the more interesting QCD pa- 
rameters estimated in Section |ITT| are illustrated in figures |23| and The relatively large 
uncertainty in the gluon condensate estimate reveals an anti-correlation with the mass. 
Qualitatively similar plots hold for the residue and continuum threshold. 
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FIG. 23. Scatter plot illustrating the association between the p-meson mass and the quark 



condensate a estimated in Section ED 
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FIG. 24. Scatter plot illustrating the association between the p-meson mass and the gluon 
condensate b estimated in Section ID. 
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TABLE III. QCD Parameters having significant roles in determining p-meson spectral proper- 
ties. Measures of association include Cramer's V, the Contingency Coefficient C, the Symmetrical 
Uncertainty Coefficient U, and the linear correlation coefficient r. 



Spectral Property 


Parameter 


V 


C 


U 


r 


p-meson Mass Mp 


im) 


0.34 


0.82 


0.23 


+0.78 






0.20 


0.63 


0.07 


-0.27 




Kp 


0.19 


0.62 


0.08 


+0.37 




as /it 


0.18 


0.61 


0.07 


+0.38 


Threshold Wp of (|5.1|) 


m 


0.36 


0.83 


0.24 


+0.80 




as/ir 


0.21 


0.66 


0.09 


+0.42 






0.19 


0.62 


0.08 


+0.37 


Pole Residue of dp) 


im) 


0.34 


0.82 


0.23 


+0.78 






0.21 


0.66 


0.09 


+0.46 




Kp 


0.19 


0.62 


0.08 


+0.36 



C. Nucleon Associations 



Table |IV| identifies the QCD parameters having significant correlations with the spectral 
properties of the nucleon obtained from an analysis of ( |6.1b|) at j3 = —1.2. Here the gluon 
condensate plays the dominant role in determining the mass and continuum threshold. The 
mixed condensate and finally the quark condensate play lesser roles. The predominant role 
of the quark condensate is to determine the residue of the nucleon pole. Fortunately, the 
quark condensate and nucleon mass are at least positively correlated. 

However, these correlations are dependent on the choice of interpolating field. Table 
identifies the QCD parameters having significant correlations with the spectral properties 
of the nucleon obtained from the analysis of ( [7.6b| ) and (|7.6c|) at /5 = 0. This time the 
quark condensate plays the dominant role in determining the nucleon mass. This role is also 



TABLE IV. QCD Parameters having significantly correlated roles in determining nucleon spec- 
tral properties via an analysis of ( 6.1b| ) at /? = —1.2. Measures of association include Cramer's V, 
the Contingency Coefficient C, the Symmetrical Uncertainty Coefficient U, and the linear correla- 
tion coefficient r. 



Spectral Property 


Parameter 


V 


C 


U 


r 


Nucleon Mass M^r 


ilfG') 


0.30 


0.77 


0.20 


+0.76 




{qga- G q) 


0.22 


0.66 


0.09 


-0.47 




m 


0.23 


0.68 


0.09 


+0.46 


Threshold W2 of (6.1b|) 


ilfG') 


0.31 


0.77 


0.22 


+0.79 




{qga- G q) 


0.22 


0.64 


0.09 


-0.44 




im) 


0.22 


0.65 


0.09 


+0.43 


Pole Residue of (p.lbj) 


m 


0.30 


0.78 


0.17 


+0.69 




i^G') 


0.24 


0.70 


0.15 


+0.63 




{qga- G q) 


0.16 


0.56 


0.07 


-0.29 
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TABLE V. QCD Parameters having significantly correlated roles in determining nucleon spec- 
tral properties via a simultaneous analysis of ( 7.6b| ) and ( 7.6c ) at /? = 0. Measures of association 
include Cramer's V, the Contingency Coefficient C, the Symmetrical Uncertainty Coefficient U, 
and the linear correlation coefficient r. 



Spectral Property 




Parameter 


T/ 
V 




u 


r 


Nucleon Mass M^r 


{qq) 


0.31 


0.79 


0.18 


-0.68 






{qga- Gq) 


0.26 


0.73 


0.17 


+0.67 








0.16 


0.54 


0.07 


+0.21 


Threshold of 


im) 


0.32 


0.80 


0.17 


-0.66 


(7.6b|) and (7.6c|) 


{qga-Gq) 


0.26 


0.71 


0.15 


+0.64 






(ifG') 


0.16 


0.55 


0.08 


+0.52 


Pole Residue \o A3/2 


(qga-Gq) 


0.34 


0.82 


0.25 


+0.82 


of (f7.6b 


) and (|7.6cD 




0.17 


0.56 


0.08 


-0.37 






i^G') 


0.16 


0.56 


0.07 


+0.33 



shared with the mixed condensate (qga-Gq). Perhaps the most interesting result is that 
the magnitude of the quark condensate is anticorrelated with the nucleon mass, in complete 



contradiction with the "loffe formula". Figures 25, 26 and 2i illustrate the correlations 



between the nucleon mass and the quark, gluon, and mixed condensates respectively. 

What this indicates is that the intimate relationship between the quark condensate and 
the nucleon mass suggested in the "loffe formula" is invalid. Considering the assumptions 
made in deriving the "loffe formula" one might expect the quark condensate value to be tied 
more tightly to the continuum model than the ground state pole, and this is borne out in 
Table |V|. There the quark condensate plays the dominant role in determining the continuum 
model threshold. 

Of course, all the condensate values themselves are correlated by the theory of QCD, 
and this correlation has not been included in this analysis. It may very well be that the 
nucleon mass is positively correlated with the quark condensate when the higher dimension 
condensates are also correlated with the quark condensate. The important point to recognize 
is that the correlation does not come about through some direct connection such as the "loffe 
formula" . 

Moreover, we are extracting nucleon properties in a Borel regime in which these leading 
order condensates just happen to be the dominant ones. One could also extract the mass 
from much deeper in the non-perturbative regime, perhaps by a lattice QCD calculation. 
There, it is much higher dimension operators which account for the same nucleon mass [§. 
Therefore, it is really impossible to identify a particular condensate for giving rise to ground 
state hadron properties. 

The only true statement that may be made about the relationship between condensates 
and spectral properties is what we have known all along: The condensates reflecting chiral 
symmetry breaking play a significant role in determining hadron properties. 
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FIG. 25. The association between the nucleon mass and the quark condensate a estimated in 
Section JTJ The nucleon mass is extracted from the simultaneous analysis of sum rules ( |7.6bD and 



(7.6c) at /3 = 0. 
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FIG. 26. The association between the nucleon mass and the gluon condensate b estimated in 
Section [III[ The nucleon mass is extracted from the simultaneous analysis of sum rules ( |7.6bD and 



(fTiecD at /3 = 0. 
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FIG. 27. The association between the nucleon mass and the mixed condensate parameter ttiq 
estimated in Section |ll|. The nucleon mass is extracted from the simultaneous analysis of sum 
rules (|g) and ( p^ ) at /? = 0. 



D. Phenomenological Correlations 



Correlations among tlie plienomenological fit parameters themselves are qualitatively 
different from the correlations with the QCD parameters. Figures ^ and 30 display the 
close correspondence between all three of the fit parameters, Wp, /p, and Mp obtained from 
1000 fits of (|5T| ). This close association between the p- meson spectral parameters suggests 
an underlying constraint on the fit parameters. Indeed, local duality necessarily leads to 
correlations among the spectral parameters. To examine the correlations of local duality 
in an analytical but qualitative fashion, consider the finite energy sum rules (FESR) 
obtained from 

so 1 fSO 

(f){s) pv{s) ds = - (f){s)lmll{s + ie)ds . (8.14) 



TT 

Here, pv{s) is given by ( p.l2b| ) and n(g^) is the QCD evaluation of ( |2.1j ) at the structure 
dfiu — P^iPu/p"^- The weighting function 0(s) for FESR is taken to be 1, s, s^, . . .. These 
functions weight excited state contributions in favor of ground state contributions and, as 
such, FESR are inappropriate for a quantitative determination of ground state properties. 
However, the first moment and, to a lesser extent, the second moment of the FESRs can 
provide some general insight into the relationships among the spectral parameters. Taking 
0(s) = 1 and 0(s) = s yields the following two relationships 

87rVp-So = 0, (8.15a) 
Svr^ K'§ = (^^) + y i^G') ^ . (8.15b) 
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FIG. 28. The association between the p-meson continuum threshold and the decay constant 
obtained from 1000 fits of (O). 



In a qualitative sense, the contributions from the OPE on the right-hand side of (|8.15b|) are 
negligible. In this case ( 8.15| ) implies 



p ' 



(8.16a) 
(8.16b) 
(8.16c) 



These relationships, capturing the features of local duality in a qualitative sense, are inde- 
pendent of the precise values of the QCD parameters. The relationships remain prominent 
in the more quantitative Borel sum rule analysis as indicated in figures |2^, ^ and Of 
course, the Borel analysis is done independent of these qualitative constraints. The first 
relation, the most reliable FESR relation, is satisfied remarkably well. 

An even more striking correlation is illustrated in figure |3l| for the nucleon sum rule 
of ( H.lb ) with the optimal (3 = —1.2. The association between the continuum threshold 
and the nucleon mass is nearly one-to-one despite the underlying scatter of QCD parameter 
values. The explanation for this phenomena is that there is a single term in (|6.1b|) which 
is capable of resolving the pole from the continuum. The two leading terms of the OPE 
are used in deriving the continuum model. Hence if the third term is vanishingly small, the 
best fit is most certainly pure continuum as the sum rule is satisfied as a perfect 



Laplace transform. The few points lying above the curve in figure ^ where the mass equals 
the threshold are cases where the third term in the sum rule has failed to resolve the pole 
from the continuum. Figure ^ displays the corresponding association for the mass and the 
residue of the pole. 
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FIG. 29. The association between the p-meson mass and the continuum threshold obtained 



from 1000 fits of (5.1 
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FIG. 30. The association between the /)- meson mass and the decay constant obtained from 
1000 fits of Q. 
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FIG. 31. The association between the nucleon mass and the continuum threshold obtained 
from 1000 fits of ( |6.1bD for the optimal /? = —1.2. The dashed line indicates the equivalence of the 
mass and threshold. 
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FIG. 32. The association between the nucleon mass and the residue of the pole obtained from 
1000 fits of ( |6lb| ) for the optimal P = -1.2. 
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In the p-meson case there are in fact two terms independent of the continuum model in 
( ^.1|) . Both terms act together to resolve the pole and continuum and hence give rise to 
some scatter in the mass-threshold association. Scatter in the residue of the pole reflects 
the fact that all terms affect the strength of the pole. 

The strong associations among nucleon spectral parameters reflect the underlying con- 
straints of local duality. Relationships analogous to the FESR for the p meson may be 
derived for the nucleon spectral parameters However, the larger continuum model con- 
tributions to the nucleon correlator make the use of FESR less desirable. The approximate 
nature of the continuum model may be too crude for a weighting function that fails to 
suppress excited state strength in the spectral density. The resultant FESR relations are 
relatively crude and do not capture the essence of local duality as well as in the p-meson 
case. 

A very different picture develops for the correlations among the nucleon fit parameters 
obtained from the analysis of ( |6.1b| ) for the unfavorable value of /? = +0.8. Figure ^ displays 
the correlation between the mass and the continuum threshold, W2- The corresponding 
correlation between the mass and the residue is illustrated in figure |3^ . 

The nature of the three regimes discussed earlier in section |VIF| is resolved here. The 
near-pure continuum fits are revealed in the lower-left corner of the mass-threshold scatter 
plot of figure 33. The collision of the continuum threshold and the nucleon mass gives rise 
to the intermediate regime where the gap between ground and excited states has not been 
resolved. The smaller residue compensating for strength from the continuum model in these 
fits is illustrated in figure 0. The third pole-plus-continuum regime displays the same near 
one-to-one correspondence between the mass and the threshold as seen in figure ^ for the 
optimal (3 = —1.2. 

The correlations among the spectral parameters are interesting because they reveal the 
manner in which the sum rules work. Correlations among the fit parameters from the 
simultaneous analysis of ( [7.6 b| ) and ( |7.6cD at /5 = are illustrated in figures and 
These distributions are similar to those for the celebrated p-meson sum rule, and reflect the 
presence of two terms separating pole from continuum. 



IX. ON THE NECESSITY OF DIRECT INSTANTONS 

Arguments favoring the need for direct instanton contributions to the QCD side of QCD- 
SRs |]7(]| , p(]| are based on the assumption that the average instanton size [|7^], pc = 0.33 fm 
<C 1/p = 0.39 fm, the separation scale of the OPE. Even though p^ ^ 1/fi, proponents of 
direct instanton contributions still argue that the effects of instantons are not sufficiently 
accounted for in QCD condensates alone. 

The analysis of direct instantons in nucleon sum rules are based on many additional 
assumptions. Here we focus on the most recent investigation of Ref. [Q. There, correc- 



tions to the identity operator are neglected. As discussed in Section these leading order 
corrections are 50% and suggest that as corrections are completely out of control. Moreover 
the possibility of a significant dimension-two power correction arising from a summation of 
the perturbative series, uncertainties surrounding the factorization of the four-quark oper- 
ator, coupled with an absence of a valid Borel regime for the sum rule of ( |6.1a| ) make any 
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FIG. 33. The correlation between the nucleon mass and the continuum threshold, 11)2, obtained 
from the analysis of ( 6.1b| ) for the unfavorable value of /? = +0.8. The dashed line indicates the 
equivalence of the mass and threshold. 
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FIG. 34. The correlation between the nucleon mass and the residue of the pole, Xo, obtained 
from the analysis of (|6.1b| ) for the unfavorable value of /3 = +0.8. 
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FIG. 35. The correlation between the nucleon mass and the continuum threshold, w^, obtained 
from the simultaneous analysis of sum rules (|7.6bD and (|7^ at /? = 0. 
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FIG. 36. The correlation between the nucleon mass and the residue of the pole, Xo^s/2i ob- 
tained from the simultaneous analysis of sum rules (|7^) and (|7^ at /? = 0. 
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conclusions dubious. 

In addition, the OPE is truncated at dimension six. As we saw in Section ^ the 
dimension six, seven, and eight operators are crucial to separating the ground state pole 



from the continuum. Without the dimension seven term, the best fit of (|6.1b|) is pure 
continuum. The fit is perfect when W2 —>■ and A ^ 0. Hence, the nucleon sum rule at 
the structure 1 in Ref. ||60| (equation 12) is unstable by design. Any further discussion of 
stabilizing the sum rule by adding instanton contributions is uninteresting. Finally, issues 
surrounding the problems of double counting were handled crudely by simply eliminating 
OPE contributions beyond dimension five. Hence one should not be prejudiced by existing 
analyses suggesting nucleon sum rules require direct instanton contributions. The issue is 
really not settled. 

For the spin-1/2 interpolator of ( |2.4a| ), direct instanton contributions to the sum rule of 
( |6.1b|) are governed by the factor fST 



C2 = (l3 + 10 /? + 13) /8 . (9.1) 

The direct instanton contributions were argued to be of vital importance in Ref. 0] in the 
vicinity of (3 ~ —1. 



In the lattice QCD investigation of ||T3[ it was discovered that the overlap of the interpo- 
lator X2 with the nucleon ground state is one hundred times smaller than that for xi- While 
the nucleon mass is certainly independent of the interpolating field mixing parameter jS, one 
also has the residue of the pole, Xq, independent of /3 at the 1% level. This anticipated 
independence of both ground state properties provides new opportunities for evaluating the 
necessity of direct instanton effects. 



As a first test of sum rule consistency, we utilize the sum rule of ( |6.1b| ) to extract the 



mass as a function of the interpolating field mixing parameter p. As /? — —1.5 the terms 
contributing to the continuum model become very small and additional poles must be added 
to the phenomenological side of the sum rule to account for strength in the correlator lying 
above the ground state mass. As a result we restrict (3 > —1.3. Valid Borel regimes were 
found over the range —1.3 < (3 < —0.9. The coefficient C2 changes by 50% over this range. 
If direct instanton contributions are indeed important at (3 ~ —1, the absence of direct 
instanton contributions in this investigation might be revealed by a significant (3 dependence 
in the nucleon mass. 

Figure |3^ displays the ratio of the nucleon mass extracted at (3 to the value obtained 
at the optimal f3 = —1.2. Uncertainties are obtained from two uncorrelated samples of 200 
QCD parameter sets. As such the central value at f3 = —1.2 need not be 1. The large 
uncertainties associated with this sum rule prevent a sensitive test of the necessity of direct 
instantons. It is clear that all ratios accommodate the ratio of 1. 

A more demanding test may be performed using the sum rules arising from the consider- 
ation of the overlap of spin-1/2 and spin-3/2 interpolators. The smaller uncertainties allow 
a more interesting examination of these issues. In addition, the continuum model terms are 
independent of (3 and provide a wider range of (3 in which a valid Borel regime exists. 

Figure addresses the dependence of the nucleon mass obtained from the sum rules of 
( [7.6b|) and (|7.6c|) on the interpolating field mixing parameter f3. Here, ratios of the mass 



obtained at finite P to P = are displayed. Sum rule consistency demands that all ratios 
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FIG. 37. The dependence of the nucleon mass obtained from the sum rule of ( |6.1b ) on the 
interpolating field mixing parameter /?. Uncorrelated uncertainty estimates are obtained from two 
samples of 200 QCD parameter sets. 



agree with one, and this is satisfied in figure ^ without direct instanton contributions. The 
uncertainty estimates are obtained from two uncorrelated samples of 200 QCD parameter 
sets. We note that a valid regime does not exist for j3 < —0.8. Hence the use of the loffe 
interpolator {j3 = —1) is a particularly unfavorable choice for these sum rules. The large 
error bar at = —0.6 reflects the small valid Borel regime used in the fit. 

The corresponding dependence of on (3 is plagued by huge uncertainties. Small 
changes in the nucleon mass which appears squared in the exponential of the phenomeno- 
logical side of the sum rules leads to large fluctuations in Xq. Hence it is difficult to gain 
any information on the possible size of direct instanton contributions. However, they are 
certainly unnecessary as all ratios encompass 1. 

In order to proceed, we consider fixing the nucleon mass to the value determined at 
/3 = in an attempt to reduce the uncertainty in A©. Figure ^ displays these ratios. The 
corresponding plot for the continuum threshold is given in figure 

For (3 < the continuum threshold tends to rise and the distribution of the residue 
remains largely independent of j3 without direct instanton contributions. For (3 > the 
continuum threshold drops to account for strength that would otherwise be absorbed by the 
rising nucleon mass as illustrated in figure ESI. With the threshold near the pole position. 



strength which might otherwise rest in the pole has been suppressed as discussed in section 
[VIII D| . Hence the drop in the residue is an artifact of our fixing the nucleon mass in an 
attempt to reduce the error bar. 

In summary, the independence of the mass and residue from \(3\ < 1 is satisfied within 
uncertainties in the absence of direct instanton contributions. As the sum rules are improved 
in the future and QCD input parameters are better determined, it will be important to 
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FIG. 38. The dependence of the nucleon mass obtained from the smn rules of ( [7.6b ) and (7.6c) 
on the interpolating field mixing parameter Uncorrelated uncertainty estimates are obtained 
from two samples of 200 QCD parameter sets. 



repeat this analysis to again evaluate the necessity of direct instanton contributions or other 
refinements of the QCD-SR approach. 

At present, direct instanton contributions are not required. This is good news indeed 
as the amount of modeling required for instanton estimates is relatively unpalatable when 
compared to that for the standard QCD-SRs themselves. In fact a preferable approach 
would be to simply raise the OPE separation scale further into the perturbative regime to 
ensure that direct instanton effects are more fully accounted. However, this does not yet 
appear to be necessary. 



X. COORDINATE SPACE CORRELATOR IMPLICATIONS 

Ground state properties have been found to be independent of the interpolating field 
mixing parameter (3. On the other hand, the leading terms of the OPE for ( |6.1| ) are very 
dependent on (3. This situation has interesting ramifications for previous analyses of co- 
ordinate space point-to-point correlation functions fH^pBUT^ where there has been some 
confusion surrounding the physics represented in point-to-point correlation functions. 

In Ref. [^] some conclusions are drawn from an incorrect interpretation of the physics 
represented in two-point correlators. These authors suggest that the behavior of a ratio of 
the interacting to free correlator can reveal information on clustering in the scalar diquark 
channel. 

It should be emphasized that the phenomenological description of the two-point function 
involves the baryon mass, the continuum threshold describing the effective onset of excited 
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FIG. 39. The dependence of Xq determined by the sum rules of ( 7.6b|) and ( 7.6c| ) on the 
interpolating field mixing parameter (3. Here the nucleon mass has been fixed to the optimal value 
obtained at /3 = 0. Uncorrelated uncertainty estimates are obtained from two samples of 200 QCD 
parameter sets. 
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FIG. 40. The dependence of the continuum threshold obtained from the sum rules of ( [7.6b ) 



and (7.6c) on the interpolating field mixing parameter /3. Here the nucleon mass has been fixed to 
the value obtained at /? = 0. Uncorrelated uncertainty estimates are obtained from two samples of 
200 QCD parameter sets. 
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states and the residue of the pole which indicates the abihty of the particular interpolator 
to excite the baryon from the QCD vacuum. Information on quark configurations in the 
nucleon wave function is absent in point-to-point correlation functions. 

The coordinate-space correlator may be obtained from (|2.10b|) or (|2.11b| ) via 



The important point here is that for large x, Il{x) is dominated by the contribution from the 
ground state pole. As such n(x) is independent of \P\ < 1 at the 1% level, or alternatively 
n(a;) is independent of |/3| < 10 at the 10% level |ll. 

The coordinate space analyses focused on a ratio of Il{x)/Il^{x) where Il^{x) is the free 



coordinate-space correlator. For the correlator at the structure 7 ■ p studied in [|62| , |63| , |72| at 
f3 = —1, the generalized free correlator is given by 

n°(x) = (5 + 2/3 + 5/5^) . (10.2) 

Since the denominator is (3 dependent and the numerator is independent of (3 for large x, 
the ratio is dependent on p. The rate at which the curve turns up as the ground state 
dominates and the free correlator continues to drop off for increasing x is dependent on the 
interpolator mixing parameter (3. The rate at which the ratio increases above 1 describes the 
relative weightings of ground versus excited states generated by the interpolating field. It has 
nothing to do with the possibility of scalar diquark degrees of freedom being predominant 
in the nucleon wave function. Moreover, a comparison with a A+ correlation function for 
example has no meaning. 



XI. CONCLUSIONS 

A rigorous procedure for extracting quantities of phenomenological interest from QCD 
Sum Rules has been established. We have performed a comprehensive analysis of ground 
state p-meson and nucleon QCD sum rules. The inclusion of uncertainty estimates for 
the phenomenological parameters has allowed us to evaluate the predictive ability and self 
consistency of the QCD-SRs. When the analysis is done rigorously, QCD Sum Rules are 
predictive. 

Some have argued that fitting sum rules appears to be more of an art than a science. 
Indeed, many of the findings here contradict the conventional wisdom of both practitioners 
and skeptics alike |73] and suggests that past arguments have been based more on rhetoric 



than sound scientific arguments. 

The analysis of thousands of QCD parameter sets reveals the true predictions of the sum 
rules. For the nucleon the results may be summarized as follows. 

• The nucleon sum rule of ( |6.1aD , traditionally taken to be the most reliable sum rule, 
fails to have a Borel regime where both the truncated OPE is reasonably convergent 
and the ground state dominates the phenomenology. 
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The nucleon sum rule of ( |6.1b|) , also obtained from spin- 1/2 interpolating fields, re- 
solves the nucleon mass to ±260 MeV. This uncertainty is more than double the value 
commonly assumed by most practitioners. 



The nucleon sum rules of ( |7.6b|) and (|7.6c|) obtained from the consideration of a spin- 



3/2 isospin-1/2 interpolator projected by the generalized spin- 1/2 isospin-1/2 interpo- 
lator are much more reliable than the traditional sum rules. They have greater overlap 
with the ground state pole relative to the continuum model, have broader valid Borel 
regimes and provide phenomenological estimates of nucleon spectral properties which 
are more stable than those obtained from the conventional sum rules. 



The nucleon mass obtained from a simultaneous analysis of ( [7.6b| ) and ( [7.6c|) at the 



optimal /3 = is resolved to ±80 MeV using QCD parameter estimates from the 
literature. 

• The dependence of the spectral parameters on the interpolating field mixing param- 
eter is in accord with expectations without resort to direct instanton contributions. 
This is compelling evidence that instanton physics is adequately accounted for in the 
nonperturbative vacuum expectation values. 

The latter point is truly exciting for the QCD sum rule approach. The amount of modeling 
required for instanton estimates is relatively unpalatable when compared to that for the 
standard QCD-SRs themselves. 

The QCD-SR typically regarded as the most reliable or stable nucleon sum rule has been 
found to be invalid. A summary of additional problems with this sum rule includes: 

• The leading order corrections are the order of 50%. As such, perturbative correc- 
tions to this sum rule are very likely out of control. 

• Summation of the perturbative series may give rise to an important dimension-two 
contribution p5| , |2l| . The size of such a correction is presently unknown. 



• Contributions from the four-quark condensate are not well known and appear to be 
large. Unlike the p-meson sum rule, the four-quark condensate is not suppressed by a 
factor of Us- 

• A regime in which the OPE is reasonably convergent while the phenomenological 
description is dominated by the ground state pole is non-existent. If one proceeds 
with this sum rule, the results are very sensitive to the model for the continuum, or 
inaccuracies in the OPE lead to incorrect results as seen here. 

Most QCD Sum Rule predictions of nucleon properties are based on this sum rule. The 
results are suspect, and should be reevaluated using the new techniques introduced here. 

The "loffe formula" directly relating the nucleon mass to the local chiral condensate has 
been found to be misleading at best. While there is a strong connection between chiral sym- 
metry breaking and hadron properties, such detailed or direct relationships are impossible 
to make and it is inappropriate to be more specific on the origin of the nucleon mass. 
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The fundamental and derivative sum rules for the p meson have also been examined. Once 
again, the popular sum rule analysis has been found to be undesirable. The use of derivative 
sum rules for the determination of p-meson spectral properties is unfavorable. Increased 
continuum model contributions, poorer OPE convergence, and loss of vital information in 
the derivative sum rule results in discrepancies between the fundamental and derivative sum 
rules. 

The associations among the phenomenological fit parameters are particularly interesting 
as they reveal how the sum rules resolve the spectral properties. We have found that 
the separation of the ground state from the excited state continuum model is completely 
determined by the high dimension operators of the OPE which are not used in formulating 
the continuum model. The practice of fixing the continuum threshold to a preferred value 
has been exposed as strongly determining the outcome of the sum rule analysis. QCD 
sum rule determinations of phenomenological spectral properties are necessarily dependent 
on the treatment of high dimension operators where a factorization assumption is usually 
invoked. Moreover, careful attention to OPE convergence is crucial to extracting meaningful 
results. These issues are somewhat unsettling, as we are working with expansions which are 
asymptotic. 

The emphasis here has been on exploring the QCD parameter space via Monte Carlo. 
Further refinements of the parameter space might be made by adjusting the valid Borel 
regime for each QCD parameter set such that the 10%-50% criteria are satisfied. This was 
not done here in order to explore other criteria and reduce the sensitivity to a reasonable 
yet somewhat arbitrary criteria. In addition, the Holder inequalities might be applied as a 
further check of the validity and consistency of the sum rulesQ |[74|| . 

We are now in an excellent position to attack phenomena which is of current interest 
to the nuclear physics community. The techniques are currently being applied to an in- 
medium study of vector meson properties 0. There the debates over the behavior of the 
vector-meson masses and the sum rules to be used in extracting vector meson properties in 
nuclear matter are resolved. Using standard in-medium condensate estimates, the analysis 
leaves no doubt that vector-meson masses decrease with increasing density. 

These techniques are also being applied to in-medium properties of the nucleon 0, by 
generalizing the vacuum sum rules analyzed here to finite nuclear matter density. The 
Monte-Carlo techniques are resolving the predictive ability of the in-medium sum rules. 

Another investigation utilizing these Monte-Carlo techniques studies the QCD-SR deter- 
mination of the p — uj meson mixing parameter A |]^. There, correlations between the ratio 
of u and d quark masses are particularly interesting. 

The relatively large uncertainty in the residue of the pole does not bode well for QCD 
Sum Rule analysis beyond two-point functions. For example, the residue of the pole appears 



as a pre-factor to any nucleon matrix element such as qa or magnetic form factors |T5|]. It 



will not be surprising to find uncertainties in these observables as large as the 35% relative 



We note that the shaded regions satisfying the Holder inequalities in Ref. [74| are determined 
using both fundamental and derivative sum rules. The regions are altered when determined by the 
fundamental sum rule alone. 
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error in the squared residue itself. As the uncertainties in the additional vacuum matrix 
elements relevant to the three-point functions are folded in the uncertainties are likely to 
exceed 50% and may in fact approach 100%. Analyses of these issues are currently underway 

ZI,lZl- 



When compared to uncertainty estimates previously published in the literature, our 
uncertainty estimates are somewhat larger. This is not due to any shortcoming of the 
approach presented here. Instead, previous authors have failed to provide a reliable method 
for determining realistic uncertainties. Often, uncertainty estimates have been conjectured 
without any supporting quantitative analysis. 

This in-depth examination of QCD sum rule self consistency paints a favorable picture 
for further quantitative refinements of the approach. Of course such progress will be very 
difficult, as there are many aspects of the approach which must be refined for any particular 
improvement to be relevant. The techniques introduced here, such as maintaining inde- 
pendence of ground state properties from the interpolating field mixing parameter, might 
be used to further narrow the uncertainties of the approach. It is extremely important to 
refine the QCD-SR approach such that uncertainties in the QCD input parameters may be 
reduced while maintaining x^/^bf — 1- Research in this direction is required before mean- 
ingful QCD-SR-based predictions for nucleon matrix elements and moments of structure 
functions are possible. 
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